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ABSTRACT 

A  direct  design  method  for  solving  the  problem  of 
robustness  to  cross-ccupling  perturbations  in  ma ltivariable 
control    systems     is    presented.  The   method      uses   numerical 

optimization  procedures  to  manipulate  the  system  feedback 
gains  as  direct  design  variables.  The  manipulation  is  accom- 
plished in  a  manner  that  produces  desired  performance  by 
pole  placement  and  robustness  by  modification  of  the  minimum 
singular  values    of    the  system   return   difference    matrix. 

Channels  affected  by  cross-coupling  perturbation  may  be 
recognized  by  the  character  of  their  transfer  function  plot. 
The  nechanism  used  by  the  pole  placement  and  robustness 
routine  in  obtaining  a  robust"  design  is  evideat  from  the 
gain  changes  associated  with  the  transfer  function  diagram 
and  the  zero  shifts  shown  on  pole-zero  plots.  The  pole 
placement  and  robustness  routine  uses  gain  equalization  and 
zero  assignment  to  modify  the  .  characteristics  of  the  system 
in  the  areas  of  low  singular  values,  producing  a  robust 
design. 

A  modification  of  the  pole  placement  and  robustness 
routine  that  may  be  applied  to  the  design  of  robust 
observers  is  also  presented.  Using  feedback  and  filter  gains 
as  direct  design  variables  a  practical  design  procedure  for 
robustness    recovery    in  observer  based   systems  is   obtained. 
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I.    IHTEODOCTION 

With  the  rising  interest  in  muitivariable  coatrol  theory 
brought  on  by  increasingly  complex  systems  the  need  has 
arisen  to  develop  design  methods  that  will  allow  the 
designer  tc  specify  system  performance  while  at  the  same 
time  ensuring  relatively  high  stability  margins  or  robust- 
ness. In  the  single-input  single-output  (SISOi  case  the 
designer  has  had  the  tools  to  do  these  tradeoffs  in  the  form 
of  Nyguist,  Bode  and  root  locus  plots.  In  the  multi-input 
multi-output  (MIMO)  case  the  classical  methods  are  not 
totally   appropriate. 

With  the  increased  interest  in  MIMO  systeas  numerous 
methods  of  design  have  been  employed  to  obtain  suitable 
system  performance  and  robustness  with  varying  degrees  of 
success.  One  primary  method  of  design  is  to  keep  the  plant 
as  decoupled  as  possible  throughout  the  design  so  that  each 
individual  element  may  be  controlled  independently  and 
designed  essentially  as  a  single  loop  system.  Rosenbrock 
[Eef.  1]  has  developed  a  procedure  where  the  multiloop 
system  is  modified  into  a  system  that  has  diagonal 
elements  that  are  much  larger  than  any  off-diagonal 
elements.  This  diagonally  dominant  system  is  then  in  a  form 
where  conventional  Nyguist  type  techniques  can  be  employed 
in  the  analysis.  A  third  common  MIMO  design  method  is  that 
of  the  linear  Quadratic  {LQ)  method.  This  method  uses  a 
guadratic  cost  functional  and  optimization  principles  to 
allow  the  designer  tc  design  for  various  performance  levels 
by  adjusting  the  matrix  weighting  terms  used  in  the  cost 
function.  The   major      difficulty      with      all  of      the      above 

methods  is  that  they  are  not  necessarily  robust .  This  is 
especially   true    for    cross-coupling   terms    between    loops. 
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The  primary  achievement  of  this  thesis  has  been  the 
incorporation  of  the  time  domain  pole  placement  design 
procedure  with  a  method  of  using  return  difference  matrix 
singular  values  to  improve  the  robustness  of  the  design.  Ihe 
technique,  which  utilizes  a  modern  optimization  routine,  can 
significantly  assist  the  designer  in  obtaining  robustness  in 
the  face  of  cross-coupling  perturbations.  It  has  also  been 
shown  that  the  cross-coupling  perturbation  problem  can  be 
detected  ty  using  classical  open-loop  Bode  diagrams  as  well 
as  modern  control  analysis.  The  pole  placement  and  robust- 
ness design  routine  developed  for  this  thesis  his  been  used 
on  several  problems  discussed  in  recent  literature.  In  these 
studies  the  pole  placement  and  robustness  design  code  has 
proven  capable  of  meeting  the  desired  goals  cf  pole 
placement  and  robustness  and  also  brought  to  light  seme 
interesting  aspects  of  the  cross-coupling  perturbation 
problem.  A  slightly  modified  pole  placement  and  robustness 
routine  has  proven  effective  in  the  design  of  robust 
observers. 

The  remainder  of  the  thesis  will  present  background 
material  on  SISO  systems  in  Chapter  Two  and  on  MIMO  systems 
in  Chapter  Three.  Optimization  will  be  discussed  in  Chapter 
Pour  along  with  a  discussion  of  the  Automated  Design 
Synthesis  (ADS)  program  used  as  the  optimizer  routine  for 
the  pc]e  placement  and  robustness  design  procedure  developed 
in  this  thesis.  The  thesis  methodology  will  be  discussed  ana 
outlined  in  Chapter  Five.  This  will  be  followed  by  chapters 
discussing  results  from  selected  problems.  Conclusions  will 
be   presented  in   the    final  chapter. 
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II.    SINGLE-INPUT    SINGLE-OOTPOT    SYSTEMS 

The  purpose  of  this  thesis  research  has  been  to  develop 
a  method  of  obtaining  a  robust  multivariable  control  system 
design.  A  brief  review  of  the  concept  of  robustness  and 
stability  in  the  framework  of  a  conventional  SISO  system 
will  be  dene  before  pursuing  the  concepts  in  a  more 
complicated  fashion  in  the  following  chapters.  A  simple 
interpretation  of  robustness  is  the  ability  of  the  system  to 
tolerate  design  perturbations.  These  perturbations  could  be 
in  the  form  of  actuator  failures,  plant  parameter  uncer- 
tainty, unmodeled  dynamics  or  nonlinear  terms,  or  any  one  of 
many  other  perturbations  to  the  nominal  design  of  the 
system. 

The  primary  reason  for  feedback  systems  is  the  control 
cf   uncertainty      within  the    system.  By   appropriate      use  of 

feedback,  properties  that  would  lead  to  an  unstable  system 
may  be  controlled.  When  stability  and  robustness  aspects 
are  considered  for  a  SISO  system,  frequency  domain  design 
concepts,  using  either  Nyguist  or  Bode  plots,  are  normally 
used.  Robustness  in  SISO  systems  is  formulated  naturally  by 
the  concept  of  gain  and  phase  margins,  both  of  which  are 
readily   available  on   the   Nyguist    or   Bode   diagram. 

In  figure  2.1a  nominal  feedback  system  can  be  seen  with 
a  perturbation  element  O^  (s)  placed  in  series  with  the 
nominal  system.  When  o£  =  1  the  system  is  nominal  and  stable, 
lo  determine  the  positive  phase  margin  of  the  system  the 
value  of  o£(s)=c((jw)  =  e  will  be  changed  sy  varying^ 
until  the  system  just  becomes  unstable.  This  value  of  </  will 
then  be  the  system  phase  margin.  The  negative  phase  margin 
can  be  computed  in  the  same  way.  To  find  the  gain  margin  the 
magnitude      value   of    c{   is  increased      until      the  system      just 
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Figure   2.1        Basic   Gain   and    Phase  Margin   Determination   Model 


becomes  unstable.  This  is  the  upward  gain  margin.  A  similar 
lower    margin  is   also   defined. 

These  gain  and  chase  margins  may  not  be  adequate  meas- 
ures cf  robustness  [Ref.  2]  because  they  do  not  account  for 
simultaneous  variation  in  both  gain  and  phase.  Therefore, 
while  large  individual  gain  or  phase  changes  may  not  desta- 
bilize the  system,  small  simultaneous  changes  in  gain  and 
phase  may  destabilize  the  system.  This  is  not  a  major  diffi- 
culty in  classical  SISO  techniques  because  the  effect  can  be 
easily   detected. 

Gain  and  phase  margin  can  be  defined  in  terms  of  the 
cpen-lcop  frequency  domain  plots  in  either  the  Bode  or 
Nyguist  format.  Figure  2.  2  depicts  a  classical  Bode  plot 
showing  gain  and  phase  margin  determination  from  the  plot. 
The  Nyquist  plot  may  also  be  used  to  obtain  this  informa- 
tion. Nyguist  criterion  states  that  if  tha  open-loop 
transfer  function  G  (s)  H  (s)  has  n  poles  in  the  right  half 
plane   and   the   limit    cf  G  (s)  H  (s)=constant   as   s  ->-©    then    for  a 
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Figure  2.2  Classical  Bode  Plot- 
stable  system  the  locus  of  G(s)H(s)  will  encircle  the  -1+jO 
point  n  times  in  the  counterclockwise  direction  as  s  varies 
along  the  Nyguist  contour.  If  there  are  no  poles  in  the 
right  half  s  plane  then  the  locus  will  not  encircle  the  -1  + 
jO  point.  The  diagram  in  figure  2.3  illustrates  a  nominally 
stable  system.  The  gain  and  phase  margin  may  be  determined 
directly   from    the  diagram. 

Any  change  in  the  loop  transfer  function,  provided  the 
order  of  G(s)H(s)  does  not  change,  that  changes  the  number 
of  times  the  locus  of  G(s)H(s)  encircles  the  (-1,0)  point  in 
the  Nyguist  plot  causes  the  system  to  become  unstable.  This 
leads    to      the  conclusion      that   the     minimum   distance      of    the 
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Reg-(jcJ) 


Figure  2.3    Nyguist  Plot  of  Stable  System, 

locus  of  G(s)H(s)  to  the  (-1,0)  point  is  a  measure  of  the 
system  stability.  Ihis  distance  concept  cirries  over 
directly  to  the  MIMC  system  as  will  be  shown  in  the  next 
chapter.  Examples  cf  a  multiplicative  perturbation  and  an 
additive  perturbation  illustrate  this  idea.  Figure  2.4  is  an 
additively  perturbed  system.  Figure  2.5  shows  the  Nyguist 
plot  for  this  system.  Assuming  that  the  plant  is  itself 
stable  and  the  perturbations  are  also  stable  the  diagram  may 
then  be  used  to  determine  how  near  the  system  is  to  insta- 
bility for  the  given  perturbation. 
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Figure    2-4        Additively    Perturbed   System. 

Since  the  system  is  stable  the  (-1,0)  point  is  encircled 
the  correct  number  of  times  by  the  nominal  plant.  If  the 
locus  of  g(jcO)  in  the  diagram  is  warped  until  it  passes 
beyond  the  (-1,0)  point  then  clearly  the  number  of  encircle- 
ments of  this  point  will  change  and  the  system  will  become 
unstable,  assuming  the  order  of  the  plant  is  not  changed  by 
the  perturbation.  To  keep  the  locus  of  points  from  moving 
beyond    the    (-1,0)    point  equation    2.1      must    hold. 
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Figure  2.5   Additive  Nyguist  Plot. 


l  ^g(jw)  i<l   1  ♦  g(joi)  i 


(2.1) 


This  condition  is  illustrated  in  figure  2.6.  Tha  right-hand 
side  of  equation  2.1  is  just  the  magnitude  of  the  return 
difference  transfer  function  of  the  nominal  system.  The 
multiplicative  case  is  depicted  in  figure  2.7  with  its  asso- 
ciated Nyguist  plot  in  figure  2.8.  The  requirement  for 
stability  is  similar  to  the  additive  case  and  may  be  stated 
in   eguation   2.2 


I  ^g(jw)  I  <l      T  ♦    (g(jaj)-Ml 


(2.2) 
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Figure   2.6        Nyguist  for    Inequality  Additive  Condition. 

The    above      arguments   will  be  applied     again  in  Chapter      3   to 
develop    irultivariable  stability   and   robustness   properties. 

The  linear  quadratic  design  has  been  the  primary  method 
employed  in  modern  control  design  practice.  In  this  method 
an  optimal  state  feedback  control  law  is  developed  to  find  a 
set  of  feedback  gains  that  optimizes  a  chosen  performance 
index.  The  performance  index  for  the  steady-state  case  is 
given   in  equation   2.3 
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Figure  2.7   Multiplicative  System- 


Pi  =  /  (xfQ  x  ♦  u'B  u)dt 


(2.3) 


where  the  x 'x  term  and  the  u'u  term  form  quadratics.  The  Q 
and  E  matrices  are  chosen  by  the  designer  to  provide  the 
test  compromise  between  the  minimum  error  of  the  system  and 
the  minimum  energy  needed  to  control  the  system.  The  LQ 
method  is  based  on  the  use  of  closed-loop  state  variable 
feedback  fcr  the  control  of  the  system.  In  the  MIMO  problem 
IQ  methods  have  been  used  extensively  because  of  their  guar- 
anteed stability  margins  with  diagonal  weightiag  matrices. 
For  diagonal  weighting  matrices  the  LQ  method  yields  a 
guaranteed  phase  margin  of  60  degrees  and  -6  db  to  infinite 
gain    aargin. 
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Figure    2.8        Nyguist  Plot   for  Multiplicative    System. 

For  SISO  linear  quadratic  optimal  regulators  these 
stability  margins  can  be  developed  from  inequality  2.4 
[Bef.    3] 


|  1    ♦   f  ■  (jcJI   -   a)  -ih|  >   1 


(2.4) 
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Writing    the   Kalman    eguation    as   2.5 


|1    ♦   f«  (sl-a)-ibl  z    =    1   +     (1/,)  I    G^Cs)! 


(2.5) 


cne   has   further 


|     1    +    G(s)H(s)  |2         =    1    +     (V^)  iGK(s)|2 


(2.6) 


and  for  all  s  =  j  to  and  0<  co  <»  the  function  0/p)GK{s)  is 
greater  than  zero,  therefore  the  Kalman  inequality  is  shown 
to   be   2.7. 


|     1    +    G(s)H(s)  |    >     1 


(2.7) 


R(S)+^-n  u(S) 


■* 


Y(S) 


Figure  2.9        Typical    Eguivalent   Feedback  Systei, 
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The  graphical  interpretation  of  this  result  is  that  the 
polar  plot  cf  G(s)H(£)  must  remain  outside  the  unit  circle 
centered   at   the      -1    +  jO   point   for      all    frequencies.      Figure 


UNIT    CIRCLE 


u)-U) 


Figure  2. 10    Eolar  Plot  of  an  Optimal  System. 

2-10  shows  such  a  polar  plot.  Since  the  optimal  regulator 
with  the  loop  broken  at  the  input  to  the  plant  (denoted  by 
the  x  in  figure  2.9  )  does  not  penetrate  the  unit  disk  about 
-1  +  jO  this  means  that  the  single  input  regulator  will  have 
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a  phase  margin  of  at  least  60  degrees  and  a  gain  margin 
tolerance  of  fifty  percent  gain  reduction  aa d  infinite 
upward  margin.  Further  discussion  of  this  property  may  be 
found  in  Anderson  and  Moore  [Bef.  4].  With  this  hasic 
review  of  the  concepts  of  stability  and  robustness  in  the 
classical  SISO  system  complete,  the  next  chapter  will  extend 
some  of  these  basic  concepts  to  the  MIMO  system. 
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III.    MOLT I  VARIABLE    SYSTEMS 

Linear  quadratic  design  has  performed  relatively  well  in 
aircraft  control  concepts  because  of  the  ability  tc  formu- 
late the  system  state  equation  and  to  quantify  acceptable 
performance  indices  for  the  system.  Industrial  applications 
for  LC  theory  have  been  less  successful  and  have  led  British 
researchers  to  look  into  forms  of  decoupled  design  methods. 
Cne  of  the  primary  methods  used  in  multi variable  design  is 
to  make  the  system  totally  decoupled.  This  method  allows 
each  loop  tc  be  designed  as  a  separate  entity  by  classical 
means.  Cne  of  the  primary  difficulties  with  this  method  is 
the  problem  of  finding  a  compensator  which  will  totally 
decouple  the  system.  The  method  also  suffers  from  the 
effects  cf  cross-coupled  perturbation  terms.  A  method  that 
does'  not  totally  decouple  the  system  but  makes  the  design 
problen  simplier  by  designing  a  compensator  that  only  causes 
the  diagonal  terms  of  the  transfer  matrix  to  be  dominant 
over  all  off-diagonal  terms  has  also  been  developed, 
[Ref-  5].  Classical  frequency  domain  techniques  are  then 
used  to  design  each  loop  of  the  system.  The  major  difficulty 
with  this  and  other  single  loop  design  techniques  is  their 
failure  to  account  for  cross-coupling  perturbation  terms 
that  may  interact  between  the  loops.  The  individual  loops 
may  be  highly  robust  in  these  designs  but  the  overall  system 
robustness  may  be  low  because  of  the  loop  interaction  not 
accounted  for  in  the  design.  This  is  the  precise  area  that 
the  singular  value  analysis  has  proven  so  beneficial  in 
design  .  Singular  value  concepts  may  be  applied  tc  conven- 
tionally designed  systems  to  assess  their  robustness.  for 
instance,  if  a  system  is  designed  by  LQ  methods  the  designer 
may    then   formulate      the   transfer    function   of     the    system   and 
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assess  tie  singular  values  of  the  return  difference  matrix 
of  the  system.  If  th€  minimum  singular  value  is  found  to  be 
low  at  some  critical  frequencies  the  designee  can  then 
modify  the  Q  and/or  E  matrices  chosen  in  the  LQ  performance 
index  and  recalculate  the  design.  In  this  iterative  fashion 
a  rotust   design    would  be  developed. 

A  generalization  of  the  SISO  Nyquist  theory  discussed  in 
the  previous  chapter  has  been  made  for  the  MIftO  problem. 
Ihis  generalization  leads  directly  to  the  singular  value 
concept.  The  generalization  is  expressed  in  the  form  of  the 
multivariable  Nyquist  theorem  which  requires  that  a  closed 
loop  stable  system  have  the  same  number  of  counterclockwise 
encirclements  of  the  origin  by  the  locus  of  the  det  (I  +  G  (ju>) ) 
as  the  number  of  open  loop  poles  that  are  unstable.  Ihis 
theoren   is    formally    stated    as; 

let    N[f (s)  ]   denote    the   number     of  clockwise   encirclements  of 
(-1,0)       by    the    locus   of   f(s)      as   s   traverses   the    contour   B 
of  figure  3.1   in  a      clockwise   sense.      The   closed-loop   system 
will    he    stable   if  and  only    if   for  all   R    sufficiently  large 

N[f  (s)]   =   -P 
where 

f(s)     =  -1    +   det[I  +  G(s)]   =      <#>u(s)/    ^ft(s)    -1    and 
P  =      the    number  of      closed   right-half    plaa e      zeros   of 
*Au  (S)  . 

The  application  cf  the  Nyquist  theorem  comes  through  the 
fact  that  a  multivariable  system  will  not  ba  robust  to 
modelling  errors  if  the  return  difference  matrix,  I  +  G,  is 
nearly  singular  for  some  frequency.  If  I  ♦  G  is  nearly 
singular  a  small  change  in  G  may  make  I  +  G  exactly 
singular.  This  causes  the  det  (I  +  G)  to  become  zero  and  the 
Nyquist  encirclement  count  to  change  indicating  an  unstable 
system.  It  is  possible  for  very  small  changes  in  I+G  to 
produce  large  changes  in  the  determinant  of  I+G.  The  matrix 
I+G 
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Pigure   3.1        Nyguist   D  Contour. 


10/s+a  9.99/s+a 

10/s+a  10/s+a 

has  determinant  0.1/  (s+a)2.  If  the  element  pi(  is  changed  by 
only  cne  percent  to  9.  9/s+a  the  determinant  becomes 
1.1/(s+a)2  which  is  a  significant  change  in  the  determinant 
value.  Therefore,  it  is  evident  that  det  (I  +  G|  is  net  an 
accurate  measure  of  how  near  the  return  difference  is  to 
singularity.  Researchers,        in      the        field     of      controls 

[Ref.  6],  [Ref,  7],  [Ref.  8],  [Ref.  9]  have  used  singular 
value  analysis  to  determine  how  near  the  return  difference 
matrix   is  tc  singularity. 

Since   the  number    cf   encirclements   of    the    Nyguist   diagram 
changes   as   f  (s)     passes   through  the  -1    point   or   when   det  (I+G) 
is   zero   it   is  important   to      find   how   near   the   return    differ- 
ence   matrix     I+G   is      to   being      singular.         This      nearness   to 
singularity  can   be   interpreted   as      closeness   of  the   matrix  G 
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to  the  critical  point,  -1.  A  quantity  which  can  be  used  to 
express  the  nearness  to  singularity  of  the  matrix  is  the 
minimum  matrix  singular  value  denoted  by  CT  .  Given  a  matrix 
A  the   singular    value   may   be    expressed   by   equation    3.  1 


»:.Mt 


C  (A)    =    min    (    >•    U*A)) 


(3.1) 


where     fa  (£*£)       is      the   eigenvalue   of   the      complex    conjugate 
transpose   of     A    times     A.         A      basic   MIMO      linear   system     is 
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Figure   3-2        Basic  Multi-input,    Multi-output    System. 

depicted  in  figure  3.2.  An  additive  perturbation  to  the 
plant  is  shown  in  figure  3.  3.  If  the  plant  is  stable  before 
the  perturbation  is  added  to  the  system  the  Nyguist  theorem 
will  te  satisfied  and  the  locus  of  GH  will  not  encircle  the 
-1,0  critical  point.  When  the  perturbation  is  added  to  the 
system      as      long   as      the      Nyguist      locus      is  not      forced     to 
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Figure    3-3       Additive  Perturbation. 

encircle  the  -1  point  the  system  will  remain  stable.  A 
sufficient  condition,  recalling  the  SISO  discussion  in 
chapter  2,  for  the  perturbed  Nyguist  plot  not  to  change 
encirclements  is  that  the  norm  of  the  perturbation  &G 
remain  less  than  the  norm  of  the  return  difference  matrix  as 
expressed   in  eguation  3.2. 

I  I   £>  G(jcJ)||<    1/|  |  (I+S)-M  ||  (3.2) 

<^0  >0  This  condition  will  guarantee  that  the  Locus  of  the 
det (I+G)  does  not  pass  through  the  -1  point.  If  the  1L  or 
Euclidean   norm      is    assumed    for      this   condition     the    eguation 
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3.2  may  te  expressed  in  terms  of  singular  values  as  equation 
3.3. 

(T(AG)<  <T  (I  ♦  G)  (3.3) 

This  result  states  that  as  long  as  the  maximim  singular 
value  of  the  perturbation  matrix  AG  is  below  the  minimum 
norm  value  cf  the  return  difference  matrix  the  system  will 
remain      stable.  The  problem      of     guaranteeing      robustness 

becomes  that  of  finding  the  largest  norm  of  the  perturbation 
guantity,  the  largest  singular  value,  for  which  the  smallest 
norm  or  singular  value  of  the  return  difference  matrix  will 
remain   ncn   singular. 

The    multiplicative   form    for   a      system   such  as    figure  3.4 
gives   the  similar  norm  equation   in   3.4 

11  AGJjcJ)|l<1/||  (I  +    .{G)-i)-Hl  (3.4) 

OJ    >0.      which   may    te    expressed   as; 

<T  (AG)     <    g  (    I    ♦    G-i    )  (3.5) 

Singular  value  decomposition  software  is  readily  available 
to  determine  how  near  the  matrix  I+G  or  1+  (G)~l  is  to 
singularity. 
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Figure   3,4        flultiplicative   Perturbation. 
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IV.    OPTIMIZATION 

The  purpose  of  this  chapter  will  be  to  briefly  describe 
several  of  the  currently  employed  optimization  techniques 
and  the  Automated  Design  Synthesis  (ADS)  program,  [Eef.  10], 
which  emjloys  these  techniques.  In  general,  optimization 
implies  finding  the  "best"  possible  solution  to  a  problem. 
In  actuality  the  best  solution  found  by  an  optimization 
technique  might  really  only  be  a  "better"  solution  to  the 
problem.  The  purpose  of  ADS  and  other  optimization  routines 
is  to  allow  a  rational  search  to  be  conducted  to  find  the 
best  possible  design.  The  techniques  of  numerical  optimiza- 
tion are  used  to  logically  vary  the  various  parameters  that 
affect   the   design   until   a  good   solution   is   found. 

As  an  example  of  an  unconstrained  optimization  prcblem 
consider  the  following  problem  developed  in  [Ref.  11]-  The 
problem   is   to   minimize  the    function 

F  (x)  =10x/*-20x/2xl  +  10x£2+X/2-2xl  +5  (4.1) 

F(x)  is  often  called  the  objective  function,  the  cost  func- 
tion or  the  penalty  function.  Since  there  are  no  conditions 
imposed  on  the  design  variables,  x;  and  x  z,  and  no  addi- 
tional liaits  imposed  on  the  overall  design,  the  problem  is 
considered  to  be  one  of  unconstrained  minimization.  Figure 
4.1  represents  this  problem  in  the  design  space.  From  the 
figure  it  appears  the  optimum  is  near  the  point  1,1. 
Calculus  may  be  applied  to  determine  the  optimum  exactly. 
Taking   the    derivatives 
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Figure   4.  1        Design    Space   for   Example   Problem. 


PF  (x)/3x,  =40x,  3-40x(  xz+2x,  -2=0 


(4.2) 


9f  (x)/3xr=-20x,2    +20xr=0 


(4.3) 


and   then     sclving   the  set  of     equations  it   is   found      that   x, 
=  1.0    and  x  x  =1.0. 

If  design  conditions  are  imposed  on  the  problem  then  the 
optimization  becomes  one  of  constrained  function  minimiza- 
tion. In  other  words,  while  the  minimum  of  a  function  is 
still   sought,        this    minimum      must   exist      within    the      limits 
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imposed  by  the  design  conditions  or  constraints.  The  mini mum 
of  an  unconstrained  function  will  not  necessarily  be  the 
same    minimum  for  a   constrained   function. 

As  another  example  ,  also  from  [Ref.  11],  a  design  is 
sought  which  gives  the  minimum  weight  of  a  particular 
column.    This  weight    is   expressed   as: 

where  />  is  the  unit  weight  of  the  material  and  k  the  cross 
sectional   area.     The    stress    in    the   column    is   given    by: 

CT-B/k    =   ?/Z  Dt  (4.5) 

without  going  into  detail  the  design  will  be  constrained  by 
the  allowable  stress  en  the  structure.  Other  constraints  , 
Euler  buckling  and  shell  buckling,  are  also  of  interest  to 
the  designer  of  this  column.  The  design  problem  is  then 
stated   as: 

minimize     W   =   /9/t  Dth  (4.6) 

for    the   constraints 

g  (1)     =cr  /ij  -    1<    0  (4.7) 


g(2)   =GV<ro-i<  o  (4.8) 


9(3)     =   G"/cs  -1*   0  (4.9) 

with    D>    10e-06    and    t    >10e-0  6. 

35 


g  (4)     =   t-D<    0 


(4.10) 


^         *  0*4  '  <""s     sefer   to   the      allowable   design   stress,      Euler 
buckling    stress    and    the     shell   buckling   stress  respectively.. 
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Note:  The  optimum 
doign  is  not  unique. 
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Figure  4.2   Design  Space  for  Column  Problem. 

Figure  4.2  is  a  design  space  diagram  for  this  problem.  It 
is  noted  from  the  figure  4.2  that  the  optimum  is  not  unique 
and  can  be  any  value  along  the  arc  noted  as  optimum  in  the 
figure. 
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The    general    optimization  problem   is   then    written   as: 
minimize   F(x)  objective  (4.11) 


subject    to   g  •  (x)  <    0   j=1,  m   inequality   constraint       (4.12) 

J 


hK(x)     =  0   k=1  ,1    equality  constraint         (4.13) 


x.  <    x.<   x.       i=1,n    side   constraint  (4.  14) 

ill'  *  ' 

where   X   =  col    (x  1,x2  ,.. .  ,xn  )    is   the   design   variable. 

The  methods  used  to  solve  this  problem  are  usually 
iterative.  After  the  establishment  of  an  initial  set  of 
variables  the  optimizer  will  update  this  initial  value  until 
the      optimum     values      are      found.  Again      /borrowing      from 

[Eef.  11  ]#  the  iterative  technique  may  be  demonstrated  by  a 
simple  exanple.  Figure  4.3  is  used  to  illustrate  this 
problem.      Given    the    initial    data   set   X°      the   formula 

Xi    =      X°      +   X*  Si  (4.  15) 

can  be  used  to  upgrade  this  estimate  of  X.  The  vector  S  is 
the  search  direction  for  the  iteration  and  the  scalar  quan- 
tity P<  is  the  distance  of  the  move  in  the  S  direction. 
Eeginning  at  X°  it  is  desired  to  reduce  the  objective  func- 
tion. The  search  for  values  of  X  that  reduce  the  objective 
function  is  made  in  the  5  direction  which  in  this  example  is 
the  opposite  of  the  gradient  of  the  function  at  point  X°. 
The  choice  of  S  could  be  arbitrary  as  long  as  it  reduces  the 
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Figure  4.3        Illustrative   Example   for  Iteration. 

value  of  the  objective.  If  S  is  chosen  to  be  opposite  the 
gradient  of  the  objective  then  the  search  would  be  a 
steepest  descent  search.  Once  the  direction  is  chosen  then 
the    scalar  must    be  found    to      minimize   the   objective    along 

this  direction  vector  without  crossing  a  constraint 
boundary.  The  technigue  for  finding  c<  is  to  evaluate  X,  the 
objective  and  the  constraint  functions  for  values  of  <*  using 
computer  code  and  numerical  interpolation  to  estimate  <x  . 
This  K  value  then  gives  the  minimum  value  for  F (X)  along 
this  search  vector.  Now  that  o<  for  this  direction  is  known 
the    new   equation   becomes: 


X2     =     Xl 


<**S2 


(4.16) 


and    now      X1    is    used      to   start   a      new    search   in     direction  £2 
and   compute   a      new   <^*that  continues   to   reduce      the    objective 
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function.  These  one  dimensional  searches  are  continued  until 
no  more  reduction  in  the  objective  can  be  found.  At  this 
point    no   further   design   improvement   is  possible. 

Optimization  techniques  do  not  always  lead  to  the  abso- 
lute cptimutt  when  applied  to  problems  of  practical  interest. 
The  reasons  for  this  could  be  numerical  ill-conditioning  of 
the  problem  formulation  or  simply  that  there  are  multiple 
solutions  to  the  problem.  Because  of  these  difficulties  it 
may  be  adviseable  to  choose  several  different  starting 
points  for  the  optiiizaticn  process  and  use  engineering 
judgement  as  to  the  design  most  applicable  to  the  problem 
under   analysis. 

Considering  the  unconstrained  case  first  where  the 
desire  is  tc  minimize  the  function  F  (x) ,  it  is  well  known 
that  F  (x)  will  have  a  minimum  where  the  gradient  of  F  (x)  is 
zero.    That   is: 

grad(F  (x)  )  =0  (4.  17) 

with    the   VF(x)     defined   as: 

grad    (F(x))  =     (  <2F  (i)/3x,  ,  . . .  ,  2F  (x)  /2x„)  '  (4.18) 

Figure  4.4  shows  why  this  is  a  necessary  condition  but  dees 
not  guarantee  a  global  minimum.  The  gradient  of  F  (x)  is 
zero  at  all  four  points  A,B/C,and  D.  However,  only  A  and  D 
are  minima.  A  would  be  the  global  minimum  for  the  function 
as  defined  here.  D  would  be  only  a  relative  minimum.  To 
check  that  the  zero  gradient  corresponds  to  a  minimum  the 
Hessian  matrix, i.e  the  matrix  of  second  partial  derivatives, 
can  be  examined  for  positive  definiteness.  A  positive  defi- 
nite Hessian  ensures  a  relative  minimum.  The  only  way  to 
prove  a  global  minimum  for  the  function  is  to  show  that  the 
Hessian  matrix  is  positive  definite  for  all  design  variables 
X.  A  test  that  is  seldom  possible  to  perform. 
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Figure  4.4   Plot  cf  Various  Points  With  Zero  Gradient. 

If  the  problem  of  minimization  is  a  constrained  problem 
the  situation  is  different  from  that  discussed  above.  Ihe 
objective  function  gradient  does  not  have  to  be  zero  at  the 
optimum.  Figure  4.5  illustrates  this  case.  Using  figure  4.5 
and  assuming  a  start  at  point  A  it  is  necessary  to  choose  a 
search  direction  S  which  will  reduce  the  objective  function 
while  not  violating  the  constraint  functions.  Any  direction 
that  will  reduce  the  objective  is  said  to  be  useable.  This 
is  seen  to  be  the  half  plane  sector  below  the  tangent  to 
F(x)  at  point  A.  If  a  half  plane  to  the  right  of  the 
tangent  to  the  active  constraint  at  point  A  is  considered 
then  the  feasible  sector  is  formed.  A  combination  of  the  two 
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conditions  gives  the  useable,  feasible  sector  where  the 
search  direction  S  must  be  chosen.  In  mathematical  nomencla- 
ture the  above  argument  may  be  stated  as: 

useable  direction  grad  (F  (X) )  «S<  0  (4.19) 


feasible  direction  grad  (  g  .  (x)  )  •  S<  0  for  all  j   (4.20) 

for  which  g .  (x)  =  0 . 

Feint  B  in  figure  4.5  shows  a  point  where  the  gradient 
of  the  objective  and  constraint  point  in  exactly  the  oppo- 
site direction.  At  a  point  such  as  B  the  only  search  vector 
S  that  meets  requirements  for  useability  and  feasibility  is 
tangent  to  the  constraint  boundary  and  to  a  line  of  constant 
objective  function.    This  condition  is   stated   as: 

grad(F(x))+r>igrad(gj(x)  )  +  2>^grad(hx(x)  )  =   0  (4.21) 

where  A>  >   0     and    At+W      -unrestricted. 

With  this  brief  background  in  optimization  complete  the 
Automated  Design  Synthesis  (ADS)  program  will  be  briefly 
discussed.  More  detail  on  the  ADS  routine  can  be  found  in 
[Ref.  10].  This  code  was  developed  as  a  follow-cn  tc  the 
successful  CCNMIN  code  [Ref.  12]  developed  by  Vi nderplaats. 
It  is  designed  as  a  black  box  optimizer  which  allows  the 
user  tc  choose  combinations  of  one  dimensional  search,  opti- 
mization algorithm  and  optimization  strategy.  These  will  be 
discussed  later.  For  the  user  with  specific  requirements  the 
code  may  be  tailored  by  parameter  modification  to  meet 
specific  requirements.  For  the  work  done  in  this  thesis  it 
was  assumed  the  user  of  the  code  has  no  detailed  knowledge 
of      optimization     and   will      want      to      use      the  code      in      the 
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Figure  4.5   Constrained  Optimization  Eiaaple. 

simplest  form.  As  such,  the  ADS  calls  made  from  the  main 
program  use  only  default  parameters  and  first  focward  finite 
difference  gradients.  Should  analytical  gradients  be  avail- 
able the  user  could  use  them  within  the  code  if  desired  with 
no  difficulty.  The  calls  made  by  the  user  to  the  ADS 
routines  specify  several  important  aspects  of  the  problem 
solution.  The  user  may  call  for  an  optimization  strategy  to 
be  used  in  the  routine.   This   is  not  required  and   its  use 
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depends  on  the  problem.  These  strategies  are  discussed  in 
[Eef.  10].  Two  strategies  used  most  often  in  the  analysis 
done  for  this  thesis  have  been  sequential  unconstrained 
minimization  using  quadratic  exterior  penalty  function  and 
the    augmented     Lagrange   multiplier      (AIM)       method.  Others 

available  include  sequential  linear  programming  and  sequen- 
tial   guadratic    programming. 

The  basic  optimizer  is  also  chosen  by  the  user  from  two 
unconstrained  and  three  constrained  optimization  algorithms. 
The  unconstrained  algorithms  are  Fletcher-Reeves  conjugate 
directions,  Davidon-Fletcher-Powell  (DFPJ  variable  metric 
method  and  the  Broydcr.-Fletcher-Goldf arb-Shanno  (BFGS)  vari- 
able metric  method.  The  method  of  feasible  directions  and 
robust  feasible  directions  are  available  for  constrained 
minimization. 

The  user  has  available  several  types  of  one  dimensional 
searches  using  Golden  Section  or  polynomial  approximation 
techniques.  The  ADS  code  has  tailored  these  one  dimensional 
search  algorithms  for  the  unconstrained  and  constrained 
cases,  allowing  the  user  to  make  appropriate  choices  for  the 
type  of  problem  to  be  solved.  Figure  4.6  shows  the  basic 
organization  of    the    AIS   program. 

Since  it  is  not  the  purpose  of  this  chapter  to  cover  the 
subject  of  optimization  in  detail  all  the  possible  routines 
in  ADS  will  not  be  discussed.  A  few  of  the  routines  found  to 
work  well  for  the  work  completed  in  this  thesis  will  be 
briefly  described,  however.  There  are  several  methods  of 
optimizing  functions  of  one  variable  or  one  dimensional 
searches.  For  instance  a  large  number  of  points,  n,  could  be 
chosen  and  the  function  F(x)  evaluated  at  each  point.  The 
point  corresponding  to  the  smallest  value  of  F  (x)  could  then 
be  considered  the  optimum  value  of  the  one  dimensional 
search.    Ihis  method    is   hit    and   miss   and    better  methods        of 
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Figure    4.6        Organization  of   ADS   Prograa. 

locating   optimum   points  can    be  used.    The   methods    used   in   ADS 
are    Golden    Section   aid   polynomial  interpolation. 

Golden  Section  search  methods  are  easy  to  program  on  the 
digital  computer  and  do  not  require  continuous  derivatives. 
They  have  a  known  convergence  rate  and  are  reliable  for 
poorly  conditioned  problems.  The  major  drawback  of  the 
Golden  Section  routines  is  the  large  number  of  function 
evaluations  require.  The  Golden  Section  algorithm  is  simply 
illustrated  through  the  use  of  figure  4.7  Assume  that  X°  and 
X1  are  known  to  be  bounds  on  the  curve* s  minimum  value. 
Also,  the  function  values  of  F  (x)  ,  F°  and  F1  ,ace  evaluated 
and  known  at  these  points.  By  picking  two  intermediate 
points   X2   and  X3    where   X2      <   X3   and      evaluating    the   function 
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Figure    4-7       Golden   Section   Diagram. 

at  these  points  a  bound  on  the  minimum  may  ba  modified. 
Since  F2  for  this  figure  is  larger  than  F3  the  point  X2 
forms  the  new  lower  bound.  The  minimum  is  now  between  X2  and 
X1  .  If  the  function  F4  at  X*  is  determined  and  shown  to  be 
larger  than  F3  then  X4  becomes  the  new  upper  bound  on  the 
minimum.  Ey  repeating  this  procedure  the  bounds  may  be 
narrowed  to  any  desired  tolerance.  The  Golden  Section  rule 
is  applied  to  this  problem  to  reduce  the  bounds  in  the 
guickest  possible  time.  By  appropriately  picking  the  values 
of  X  * s  at  which  each  function  evaluation  is  made  an  effi- 
cient  algorithm    that    uses  the   ratio  proportion  of   the   Golden 
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Section, i.e.  X2/Xl  =  1.62803  is  developed.  By  choosing  a 
value  T  based  on  this  Golden  Section  rule  where  t  =  0.38197 
estimates   for  interior   points    X2    and   X3    can   be  made   as; 

X2    =     (1    -    £  )  XO    ♦  £-X*  (4.  22) 


X3    =    ?  xo   ♦     (1    -     t)  Xi  (4.23) 

As   each     new  bound      is   found  the      process   repeats       until   the 
accepted   level    of   convergence   is   reached. 

The  polynomial  interpolation  method  is  accomplished  by 
first  fittirg  a  polynomial  curve  to  the  points  about  where 
the  minimum  is  desired  and  then  finding  the  minimum  of  the 
polynomial  function.  For  example  if  the  function  F(X)  is 
approximated  by    a  quadratic    as: 


F    =    a  c  +a  (     X   +    a  z    X2  (4.  24) 

Then  the  value  of   X,   X    ,   where  Ff   is  zero  can  be  shewn 
to  be  : 

X  *  =  -a,  /2a  z  (4.25) 

If  a  -   is  positive  then  F  will  be  minimum.  Other  degrees  of 
polynomials  may  be  used  in  similar  fashion. 

Now  that  the  basic  one  dimensional  search  methods  have 
been  reviewed  the  next  step  is  to  examine  the  basic  opti- 
mization routines.  First,  the  unconstrained  case  will  be 
reviewed.   The  optimum  X  is  at  the  point  where 
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grad  (   F  (X*)  )=   0  (4.  26) 

Several  zero  order  methods  exist  for  the  purpose  of  finding 
the  minimum  value.  Ihese  include  random  search,  Powell's 
method  and' Box's  method.  Since  ADS  does  not  use  these 
methods    they  will  not  be  discussed. 

Automated  Design  Synthesis  makes  use  of  first  order 
methods  which  will  new  be  discussed.  The  steepest  descent 
method  is  best  known,  but  poor  in  performance.  Just  as  the 
name  implies  a  search  direction  is  chosen  opposite  the 
gradient   of   the    objective  function.    At  iteration    0 

SO   =      -   grad  (   F(X°)  )  (4.27) 

Figure   4.8      shows  this  algorithm    geometrically.  Note   that 

the  methed  simply  stair  steps  its  way  down  the  "hill"  to  the 
valley  or  minimum.  AES  uses  the  Fletcher-Reeves  modification 
to  steepest  descent.  In  this  routine  a  conjugate  direction 
is  chosen  to  improve  the  speed  of  the  search.  Figure  4.9 
shows    how   this    method  tracks   to    the    solution. 

The  variable  metric  methods  listed  earlier  are  usually 
more  powerful  than  Fletcher -Reeves  because  they  stcre  infor- 
mation that  allows  the  algorithm  to  approximate  the  inverse 
of  the  Hessian  matrix  or  second  derivative.  For  further 
discussion  cf  these  methods  see  [Ref.  11]  or  other  similar 
optimization  texts. 

AES  employs  two  direct  methods  for  constrainad  minimiza- 
tion. One  method  is  that  of  feasible  direction  and  the  ether 
is  the  methed  of  robest  feasible  direction.  Since  these  two 
methods  were  seldom  employed  within  the  work  presented  in 
this    paper   they    will   not  be    discussed. 

The  methods  chosen  to  handle  the  constrained  minimiza- 
tion prohlems  formulated  in  this  thesis  are  referred  to  as 
Sequential     Unconstrained      Minimization      Techniques       (SDMT) . 
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Figure  4-8   Steepest  Descent  Algorithn. 

SOMT  methods  are  methods  which  formulate  the  otjective 
function  and  the  constraint  functions  into  an  augmented 
objective  function  and  then  solving  the  problem  as  if  it 
were  an  unconstrained  optimization  task.  &.DS  employs 
interior  and  exterior  penalty  function  techniques  as  well  as 
an  Augmented  Lagrange  Multiplier  (ALM)  technique. 

The  exterior  penalty  function  method  is  incorporated  by 
forming  a  penalty  from  the  constraint  equations.  This 
penalty  is  of  the  form  : 


i 


P  (x)   =21  max  (0,  g  .  (x) )  2  +  «£!  (h  (x) )  2 


(4.  28) 
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Figure  4.9   Fletcher-Reeves  Conjugate  Direction  Algoritha. 

Prom  this  equation  it  is  obvious  that  P (x)  is  zero  if  all 
the  constraints  are  satisfied.  That  is,  if  g  .  (x)  <  0  and  h^ 
(x)  =  0,  then  all  conditions  of  the  penalty  function  are 
satisfied.  If  an  element  of  the  penalty  function  is  violated 
then  the  penalty  increases  as  the  square  of  the  violated 
constraint.  A  pseudc  or  augmented  objective  function  is 
formulated  where  : 


^(x,rp)=  F(x)  ♦  r^  P  (x) 


(4. 29) 


The  constant  r^  is  a  weighting  factor  for  the  penalty.  It  is 
adjusted  with  ADS  as  the  optimization  proceeds  to  allow  the 
program  to   systematically  converge  to  an  optimum  solution. 
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One  disadvantage  of  this  method  is   that  if  the  optimizer  is 

stopped   short  of   the  optimum   the   design  will   be  in   the 
infeasible  region  and  probably  not  useful. 

Criginally   developed   as   a  method   to  solve   equality 
constrained  problems  the   version  of   AIM  in  ADS  has   been 

modified  to   work   with   both   equality  and   inequality 
constraints.  The  statement  is  as  follows: 

minimize  F (x)  (4. 30) 


subject  to  h„  (x)  =  0  k=1  ,1  (4.31) 

Next  a  lagrangian  is  created  such  that 

t   . 
L(?,M  =  F(x)  +  ZAK   h<(x)  (4.32) 

Then  the  equation  is  augmented  with  an  exterLor  penalty 
function  such  that 

A(x//X/rp)=   F(X)     *^_(^uhK(x)    +   rp(h^(x))2    )  (u.  33) 

The  power  of  this  method  is  that  in  theory  precise 
constraint  matching  is  possible  whereas  in  the  exterior 
penalty  function  method  it  is  not.  The  full  details  of  the 
method  will  not  be  covered  here,  however,  the  final  form  of 
the    objective   function   will   be  included: 

A(x,A,rp)=F(X)*XOi^*i>^)*S(>li^(X)     *r     (h     (x))2     )   (4.34) 

where  S^j  =max  (g  j  (x)  ,  -Vj/2rp  j#  ana  the  update  formulas  for 
the    Lagrange  multipliers   are: 
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^?{=K   +    2rp(max(gi(x)#->j/2rf)J    j=1,m  (4.35) 


^kL  =  K>J   2rPh  *{x)      k=1'x  (4*36J 

flith  this  trief  review  of  optimization  and  the  ADS  program 
completed  the  next  chapter  will  discuss  the  program  develop- 
ment   for   this   thesis   research. 
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*•    9£2IJ3IZATION    P2SIGN    PROCEDURE 

Tc  accomplish  the  task.  of  designing  a  control  system 
with  acceptable  time  domain  performance  and  robustness 
characteristics  in  a  straight  forward  manner  a  numerical 
optimization  technique  was  chosen  as  the  method  of  iiiplemen- 
tation  of  the  design  algorithm.  Using  numerical  procedures 
to  adjust  selected  design  variables,  in  this  case  the  feed- 
lack  and/or  filter  gains,  a  desired  level  of  performance  can 
he  achieved.  This  level  of  performance  is  actually  a  combi- 
nation of  time  domain  performance  and  robustness  or 
frequency  domain  performance.  By  establishing  the  criteria 
for  the  system  performance  in  terms  of  an  optimization 
objective  and  constraint  functions  a  versatile  procedure  can 
he  developed  to  set  the  system  feedback  gains  and  affect  an 
acceptable  design  in  terms  of  performance  and  robustness. 
The  pole  placement  and  robustness  (POPLAR)  design  procedure 
uses  pole  placement  to  establish  a  designer  selected 
performance  level  and  then  a  minimum  singular  value  level  to 
establish   the   robustness. 

The  pole  placement  portion  of  the  procedure  will  be 
discussed  first.  The  pole  placement  technique  was  chosen 
because  it  was  relatively  easy  to  implement  through  a  numer- 
ical optimization  routine.  By  using  this  numerical  procedure 
it  is  also  simple  to  incorporate  robustness  into  the  proce- 
dure along  with  the  performance  requirements.  A  numerical 
technique  similar  to  one  posed  in  £Ref.  13]  was  chosen  for 
the  pole  placement  algorithm.  An  unconstrained  optimization 
routine  from  the  IBM  IMS!  library  was  used  for  this  program. 
The  routine,  a  Newton  method,  uses  adjustments  to  the  output 
feedback  gains  to  reduce  the  size  of  an  objective  function. 
This    objective   function     was   expressed  as   a      function   of   the 
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pole  location  of  the  system  and  as  the  objective  was  reduced 
the  poles  were  moved  toward  the  desired  locations.  To 
provide  more  versatility  in  the  pole  placement  algorithm  a 
method  that  can  use  constraint  functions  as  well  as  uncon- 
strained optimization  was  chosen  for  this  program.  Ihe 
designer  may  use  either  objective,  constraint  or  a  combina- 
tion of  functions  to  secure  the  desired  pole  locations.  As 
currently  iaplemented  in  the  program  the  cost  or  objective 
portion  of  the  pole  placement  procedure  is  constructed  as 
equation   5. 1 

OBJ   =21   (Afc. ->*;   )2   ♦    (Ajfc    -Xz.    )2  (5.1) 

where  ^     =    real   eigenvalue 

^i     =    imaginary  eigenvalue 

Xe    =    desired  eigenvalue   location 

^IQ  =    desired   eigenvalue   location 
The   constraint      formulation  is  a      function  tiat      must   be 
kept    negative  or   the    constraint      is   violated.      It    is   written 
as   eguation   5.2 


g(j)   ■  V  O*.-^  )2  ♦  (Xjos-^ii*2  "  r  <5-2) 

where  r  is  a  tolerance  circle  established  as  a  function  of 
pole  placement  position.  Since  the  aim  of  the  optimizer  is 
to  keep  g  negative  any  time  the  }\  function  of  the 
constraint  is  greater  than  r  the  constraint  will  become 
active,  i.e.  violated.  The  optimizer  will  then  attempt  to 
move  the  constraint  to  the  inactive  status  by  adjusting  the 
design   parameters  of    the  system. 

Consideration  of  implementation  of  the  frequency  domain 
or  robustness  portion  of  the  design  procedure  begins  with 
the  concept  of  HIHO  phase  and  gain  margins.  Several  useful 
theorems   on   singular   value    analysis   of   multiloop   systems   are 
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presented   in  [Ref.    8]-  One   of    these   theorems      relates    the 

matrix  singular  value  of  the  return  difference  function  tc  a 
parameter,  OC,  and  further  shows  that  as  long  as  the  maximum 
singular  value  of  the  perturbation  function  (L -1  -  I) 
remains  less  than  this  c<,  the  system  remains  stable.  The 
value  of  o(  is  then  related  to  gain  and  phase  margins  of  the 
MIMO  system.  The  relationship  developed  is  given  in  equa- 
tions  5. 3   and   5. 4: 

gain   margin   =    GM    =1/(1  +  *o)  (5.3) 


phase   margin    =   Ptf    =  ±cos  -»  (1-<*02/2)  (5.4) 

provided  that   equation   5.5    holds. 

?{!+£)>     *o  (5.5) 

fcr   seme     ^<0<  1 
These   phase     and   gain  margins     are   guaranteed   in      every   loop 
simultaneously. 

Universal  gain  and  phase  margin  curves,  [fief-  14],  tased 
on  the  linimum  singular  values  of  the  return  difference 
matrix   are    developed    from  equation   5.6. 


(I-i-I)  =max-/  (1- 1/k^)2  +  2/kA(1-cos^)  (5.6) 

for    all   n      with    k„  >    0.  These   curves   shown  in      figure   5. 1 

allow  the  designer  tc  pick  a  singular  value  that  corresponds 
to  a  specific  gain  and  phase  margin  for  a  given  system.  In 
addition  tc  the  universal  gain  and  phase  plot  [fief.  15] 
developes  an  optimi2er  solution  for  formulating  a  robust 
controller    using   the    CONMIN    algorithm   [fief.    12]. 
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Figure  5-1    Universal  Gain  and  Phase  Singular  falue  Plot. 

Since  the  universal  curve  in  figure  5.1  provides  a 
convenient  method  of  specifying  gain  and  phase  margins  in 
terms  of  singular  values  the  robustness  portion  of  the  pole 
placement  and  robustness  design  procedure  uses  the  minimum 
singular  value  level  of  the  return  difference  matrix  to 
determine  the  robustness.  The  minimum  singular  value  level 
is  formulated  as  an  objective  or  constraint  function  in 
equation  5.7 


J=Z.(max  (0,  (Ofc-tf"  (j»,p))  )2 


(5.7) 


The  optimization  procedure  may  be  used  to  change  feedback 
gains  until  the  minimum  singular  value  is  raised  above  this 
desired  design  level.  Although  the  same  formulation  can  be 
used  as  a  negative  constraint  function  it  has  not  been 
implemented  as  such  within  this  program.   There  are  numerous 
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ways  the  singular  value  formulation  could  be  implemented 
within  the  program  by  changes  of  the  code  if  design  require- 
ments  forced  such  changes. 

The  pcle  placement  and  robustness  design  program  is 
based  on  the  ADS  code  discussed  in  Chapter  4  to  implement 
the  design  variable  selection  procedures.  The  pole  placement 
and  robustness  program  consists  of  two  separate  programs. 
The  first  program  is  used  to  provide  designs  for  state  or 
output  feedback  problems  while  the  second  program  is  used 
for    observer  or    filter   designs. 

For  the  state  or  output  feedback  design  program  the  user 
must  input  the  plant  matrices  A„3,C  and  initial  starting 
values  for  the  feedback  matrix  F.  The  matrices  correspond  to 
the    following   linear    differential   system: 

x    =   Ax   ♦    Bu  (5.8) 


y    =   Cx  (5.9) 


u    =   -Fx_  (5.  10) 

A  feed-forward  matrix  has  not  been  considered  in  the  current 
program.  A  feed  forward  matrix  could  be  added  t3  the  proce- 
dure   if   required   for   specific   design    cases. 

As  the  design  program  is  currently  coded  the  user  may 
run  output  feedback  or  state  feedback  design  by  specifying 
the  C  matrix  as  the  diagonal (I)  matrix  for  state  feedback. 
The  program  relies  on  initial  starting  values  of  the  feed- 
back gains,  F.  As  discussed  in  Chapter  4  there  is  no 
guarantee  that  the  optimum  found  by  the  procedure  each  time 
is      the    global      optimum   or      that   the      procedure  will     always 
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converge  to  an  acceptable  solution-  The  ability  tc  select 
acceptable  starting  values  for  the  feedback  gains  will  make 
the    procedure      more    efficient      in  operation.  As    currently 

employed,  the  program  is  used  to  obtain  pole  placement  and 
robustness  for  a  given  set  of  starting  gains  and  a  selected 
optimization  routine  from  the  ADS  program.  If  the  optimizer 
is  net  ahle  to  meet  the  desired  design  goals  on  this  program 
run  two  options  are  available.  First,  change  to  a  different 
optimization  routine  from  the  list  of  available  ADS  routines 
and  rerun  the  problem.  This  was  usually  successful  in 
improving  the  design.  Second,  the  designer  uses  a  new  set  of 
starting  values  for  the  feedback  gains  and  repeats  the 
design  procedure.  Both  options  might  be  used  on  particularly 
difficult  cases. 

The  pole  placement  and  robustness  design  procedure  has 
consistently  been  able  to  find  improved  designs;  however  the 
program  does  not  always  yield  acceptable  designs.  Certain 
problems  require  changes  in  the  optimizer  routine  and  modi- 
fication in  the  initial  feedback  gain  starting  values  in 
order  to  obtain  acceptable  designs.  tJsing  the  IBM  3033  time 
share  system  the  pele  placement  and  robustness  routine 
requires  about  ten  CPU  seconds  to  work  a  second  order 
problem  and  on  the  order  of  15  to  60  CPU  seconds  to  run  a 
forth  order  problem.  The  actual  amount  of  time  varies  with 
optimization  requirements  and   time  share   utilization. 

The  observer  robustness  design  program  requires  two 
passes  of  the  ADS  program.  In  the  first  pass  the  feedback 
gains,  F,  of  the  controller  are  computed  to  ottain  the 
desired  pole  locations.  The  second  pass  of  the  ADS  routine 
is  used  to  adjust  the  observer  gains  to  recover  the  system 
robustness.  The  two  pass  procedure  was  chosen  because  it 
allows  a  smaller  number  of  design  variables  at  each  stage  of 
the  optimization  and  much  more  efficient  computer  usage. 
Figure    5.2      shows  how      the    observer      is    implemented.  This 

diagram   is    algebraicly   stated   as; 
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X     = 
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=    A    -    FC 

(5.  11) 


where  x  is  the  state,  x  is  the  estimator  variable,  F  is  the 
feedback  gain  and  K  is  the  observer  gain.  The  design  cf  the 
feedback  gains  and  the  observer  gains  are  accomplished  as 
separate  quantities  in  keeping  with   the  separation   prin- 
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Figure  5.2   Observer  Implementation. 

ciple.  In  using  tie  pole  placement  and  robustness  design 
procedure  for  the  observer  system,  initial  vala es  of  the  F 
and  K  matrices  must  te  input.  The  same  or  different  opti- 
mization techniques  from  ADS  may  be  employed. 

The  routines  contained  in  this  thesis  have  been  based  on 
the  input  additive  singular  value  level.  The  pole  placement 
and  robustness  design  algorithm  computes  input  additive, 
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output  additive,  input  multiplicative,  and  output  multipli- 
cative singular  values.  Any  of  these  singular  values  can  be 
incorporated  into  the  objective  or  constraint  formulations 
tut    have  net  been  for  this    version   of    the   program. 

In  summary,  the  pole  placement  and  robustness  design 
procedure  is  a  straight  forward  numerical  optimization 
procedure  for  the  practical  application  modern  MIMO  system 
analysis.  The  new  aspects  of  the  procedure  are  the  implemen- 
tation of  both  pole  placement  and  robustness  criteria  within 
the  same  design  program.  The  versatility  of  the  pole  place- 
ment and  robustness  design  is  obtained  by  incorporating  a 
state  of  the  art  optimizer  routine  ADS,  with  currently 
available   singular      value  computation    routines.  Using   the 

optimizer  format  fcr  the  pole  placement  and  robustness 
design  gives  the  designer  the  ability  to  modify  variables 
directly  that  affect  both  the  time  domain  or  performance  of 
the  system  and  the  frequency  domain  or  robustness  of  the 
system.  The  numerical  optimization  incorporated  into  the 
pole  placement  and  robustness  design  is  flexible  enough  to 
incorporate  the  double  pass  design  technique  using  the  sepa- 
ration of  the  feedback  and  filter  gains  to  obtain  robustness 
recovery  for  observer  based    controllers. 
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VI.  INTRO  DOCTOR Y  PROBLEM 

The  purpose  of  this  introductory  problem  is  to  review 
some  results  of  modern  multivariate  robustness  theory  using 
a  simple  problem.  This  same  simple  system  is  then  used  as  a 
test  problem  for  the  optimization  technique  developed  for 
this  thesis.  The  problem  provides  excellent  insight  into  the 
cross-coupling  problem  and  demonstrates  how  effectively  the 
pole  placement  and  icbustness  design  procedure  can  be  in 
dealing  with  the  cross  perturbation  terms. 

The  problem  chosen  for  this  introductory  analysis  comes 
from  [Ref.  9]-  Figure  6.1  is  a  diagram  of  this  tasic 
system.  In  this  problem  a  simple  plant  is  specified  by  the 
following  linear  system: 
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x2 
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y1   =  x1  (6.2) 

y2   =   x2 

A  feedback  compensation  of  the   form  of   equation  6.3   was 
assumed. 
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fihich    gives   a  closed-loop  system,    equation    6.4 


(6.3) 


x1 
x2 


-2 


0  -2 


m                   m 

«h 

XT 

* 

2 

x2 

4 

»               t 

1 


UC1 

uc2 


This      system  has     eigenvalues      at  -2      ,-2      and  is 
stable.    Using   equation   6.5 


(6.4) 


therefore 
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G=  C(sI-S)-l£ 
the  transfer  matrix  may  be  written  as  equation  6.  6 


(6.5) 


G    =  fl/s*1        b,c/s  +  ll 
[        0  1/S+1      J 


(6.6) 


ihich  gives  the  return  difference  matrix,  equation  6.7, 


I  +  G(s)    = 


s+2/s+1      r^/s+1 
0  £+2/s+  1 


(6.7) 


Figure    6.1        Basic  Hulti-input   Multi-output  System. 

The    problem    shows  the      inadequacies   of   classical   methods 
in   establishing    the    robustness   of   the   system.    A   brief   review 
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cf  the  results  will  be  presented.  Using  the  return  differ- 
ence matrix  (I+£),  the  determinant  may  be  written  as 
equation  6.  8. 


det  (I*G)  -1=2s+3/  (s  +  1)  2 


(6.8) 


Therefore,        the   multivariate      Nyguist   diagram      will    be     as 

Im 


-t 


>Re 


Figure   6.2        Hultivariable   Nyguist  for   2s+3/ (s*1)  2. 


shown   in     figure   6-2.  The   diagram      does   not     encircle   the 

(-1/0)  point  and  is  indicative  of  a  closed-loop  stable 
system.  Considered  as  a  SISO  system  one  has  a  gain  margin  of 
-1/3  to  co  and  a  phase  margin  of  ±106  degrees.  Under  this 
criteria  one  can  conclude  that  this  is  a  good  design.  This 
will    be   shewn   later   net  to    be   the   case. 

Multi-input  systeis  are  often  designed  a  loop  at  a  time. 
Applying  loop  design  to  the  system  the  transfer  function  for 
either   lcop   taken   separately   becomes   equation   6.9 
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G=1/s  +  1  (6.9) 

figure  6.4  shows  the  Nyguist  diagram  of  this  problem.  Using 
this  Nyguist  diagram  the  system  is  indicated  to  he  stable 
and  have  phase  and  gain  margins  of  GM=(-1,<»)  and 
FM=(±180).  This  analysis  does  not  show  the  true  nature  of 
the  robustness  of  the  system.  Since  the  factor  b  is  not  a 
parameter  in  either  of  the  Nyguist  curves  it  plays  no  part 
in   the   stability   determination   using   these   diagrams. 

Using  the  criteria  of  singular  values  discussed  in 
Chapter  3  a  measure  cf  the  rearness  to  instability  for  this 
problem  may  be  obtained  by  plotting  the  minimum  singular 
value   of    (I+G)  ,  (T .  (For  numerical   calculations    a      value  of 

b/2  =    50      is   assumed).  Figure    6.3    shows      the   plot      of   this 

value  vs.  frequency.  This  gives  a  minimum  singular  value  of 
about  -23  db  or  0.07  1  which  corresponds  to  a  gain  margin  of 
about  0.93  to  1.08  and  a  phase  margin  of  ±4.1  degrees.  These 
phase  and  gain  margins  are  quite  small  and  are  evident  in 
the   cress-feed    perturbation    problem. 

A  perturbed  system  as  shown  in  figure  5.5  can  be 
produced  which  will  lead  to  stability  problems  with  very 
small  values  of  perturbation.  The  closed-loop  system  if 
perturbed  by  a  small  perturbation,  5/bJZ  ,  where  b  ,t  is  a 
large  number,  will  have  as  a  characteristic  equation,  equa- 
tion   6.10 

(sI-A)=    s2    +    Us    +S  (6.  10) 

with    the      eigenvalues  of   s=-2±  V5.  There    is      one    positive 

root    in   this  solution  and  the   system   is   unstable. 

To  determine  the  nature  of  the  robustness  of  this  system 
the    return      difference   matrix   must     be   considered.  If   the 

return   difference   matrix  of    the   transfer   function     (I+G(jw)) 
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Figure    6*3        Hiniaua  Singular  Yalue  Plot,    Example   Problem. 

is  nearly  singular  at  some  freguency  co0,  then  the  multivari- 
ate systen  will  not  be  robust  with  regard  to  any  modelling 
errors  within  the  system.  This  is  because  any  small  change 
in  G  ( jw)  can  then  lake  I+G  ( joj)  singular  and  the  det  (I+Je) 
becomes  zero,  thus  changing  the  encirclements  of  the  Nyguist 
stability  pcint  and  indicating  a  system  instability.  Using 
equation   6.1  the   pole   placement   and    robustness  design   method 
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Figure  6.4       Nyguist  Diagram   of    1/s+1- 

can  fc€  demonstrated  -  The  first  step  is  to  establish  a  tase- 
line  for  the  design.  The  state  unity  feedback  model  in  this 
simple  problem  gives  pole  locations  of  -2  and  -2.  This  set 
of  eigenvalues  -  2,  -2  were  chosen  as  the  baseline  for  the 
system.  Since  a  second  order  system  requires  only  two  feed- 
back gains  to  place  the  poles,  the  diagonal  feedback  gains 
were  chosen  for  pole  placement  purposes.  The  pole  placement 
and  robustness  design  program  was  then  used  to  place  the 
poles  and  in  the  process  obtained  the  gains  required  to  do 
this  (1/1).  A  plot  cf  some  of  the  singular  value  criteria 
obtained      is  shown      in  figure   6.6.        Pole-zero     plots   of   the 
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Figure  6.5   Perturbed  System. 

closed-  loop  transfer  functions  of  the  closed-loop  transfer 
matrix  are  shown  in  figure  6.7.  In  this  case  the  pcles  and 
a  zero  are  clustered  about  the  -2  point  and  in  the  input  two 
to  output  one  channel  a  zero  is  located  at  the  -1  point  on 
the  pole-zero  diagram.  This  point  closely  corresponds  to  the 
minimum  singular  value  frequency.  The  only  significant 
observations  are  the  ability  of  the  pole  placement  and 
robustness  routine  to  place  the  poles  and  the  relatively 
poor  singular  values  indicative  of  low  robustness. 

Since  there  is  a  requirement  for  two  feedback  parameters 
to  set  the  poles  of  a  second  order  system  there  should  te  no 
additional  freedom  in  design  to  account  for  robustness. 
Case  two  was  a  run  tc  demonstrate  this  fact.  Again  allowing 
only  two  design  variables  for  the  problem  the  pole  placement 
and  robustness  desigr  program  was  run  but  with  an  objective 
function  formulated   to  adjust   the  singular  value  above  a 
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Singular  Value  Plot  for  Simple  Problem. 
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Figure   6.7       Closed-loop  Poles   for  Pole  Placement  Only   Case. 
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base  level  of  0.6.  This  level  would  correspond,  using  the 
universal  gain  margin  chart,  to  a  gain  margin  of  -4  db  to  3 
db  and  a  phase  margin  of  ±  35  degrees.  With  only  the  two 
design  variables  to  work  with  the  pole  placement  and  robust- 
ness desigr  program  was  unable  to  place  the  poles  and  adjust 
the  singular  value  level  to  the  required  value.  The  poles 
were  placed  at  -2,  -2  but  the  singular  value  minimum  was 
still    en   the     order   of  -23    db    or    0.07.  The   pole-zerc  plot 

remained      almost   unchanged.  It   is      clear   that      additional 

degrees  of  freedom  for  the  pole  placement  and  robustness 
design  program  must  be  opened  if  robustness  is  to  be 
accounted  fcr. 

Case  three  was  then  run  on  the  pole  placement  and 
robustness  program  by  adding  on  additional  degree  of 
freedom.         This      case  used      the   f  gain   as     the    additional 

design  variable.  A  good  choice  as  will  be  seen.  Allowing 
the  optimizer  routine  the  extra  freedom  to  adjust  the  addi- 
tional feedback  gain  term  an  excellent  design  was  found. 
The  singular  value  minimum  became  0.885  as  shown  in  figure 
6.8.  Using  the  universal  gain  margin  chart  this  corresponds 
to  gain  and  phase  margins  of  -6  db  -  18  db  and  ±52  degrees 
respectively.  This  is  a  considerable  improvement  ever  the 
original  design.  The  factor  that  changed  the  design  was  the 
upper  diagonal  feedback  term  which  provides  a  cancelling 
factor   for   the      cross-coupling    term  b ;-    =    50.  This   can   be 

seen   by    looking    at    the   system   matrix   equation   6.11. 


A,    =    3    -    BF 


(6.11) 


where   eguation    6.12    gives  A 


A      = 

-c 


-1-f 


'tlt'b>2  f22. 


-1-f 


Zl 


(6.  12) 
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Figure   6.8        Singular   ?alue  Plot   Case  Three. 
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If  the  upper  diagonal  term  of  the  matrix  can  be  driven  to 
zero  the  system  will  be  decoupled  into  a  diagonal  system. 
To  place  the  poles  cf  this  system  at  -2,-2  the  diagonal 
feedback  gains  must  he  one.  The  optimizer  feedback  gains 
for  this  case  were  f„  =1.00  006,  fx^=0. 99998  and,  most  impor- 
tantly, f,- =-51. 8340 1.  Using  these  values  in  the  A  matrix 
A      =      -2.00006  1.83501 

0.0  -1.99998 

The  value  in  the  upper  right  position  in  the  system  matrix 
has  been  lowered  considerably  from  the  value  of  near  50  that 
appears  in  this  position  in  the  low  singular  value  cases, 
lowering  this  system  gain  value  decreases  the  cross-coupling 
perturbation  effects  en  the  system.  This  can  be  more  graphi- 
cally demonstrated  by  figures  6.9  and  6.10.  In  figure  6.9 
the  transfer  function  shows  a  high  gain  of  approximately  35 
db  and  a  bardwidth  of  50  rad/sec.  In  figure  6.10  this  gain 
has  been  reduced  to  6.0  db  with  a  bandwidth  of  almost  2 
rad/sec.  Two  things  are  indicated  by  the  figures,  one,  the 
open-lcop  Ecde  plot  of  the  transfer  function  of  the  cross- 
coupled  channel  can  be  used  to  indicate  the  robustness 
problem  as  evidenced  by  the  high  gain  and  bandwidth  relative 
to  the  other  feedback  gains  and,  two,  the  mechanism  used  by 
the  pole  placement  and  robustness  design  procedure  to 
recover  robustness  is  to  reduce  the  relative  gain  and  asso- 
ciated bandwidth  within  the  affected  channel.  Figures  6.11 
and  6.12  which  are  for  the  input  one  to  output  one  channel 
are  approximately  the  same  as  are  the  Bode  plots  for  the 
input  two  to  output  two  channel  which  are  not  shown.  These 
figures  indicate  that  no  problem  exists  in  the  diagonal  or 
direct  coupling  terms.  The  pole-zero  diagram  for  optimiza- 
tion run  case  three  is  included  in  figure  6.13.  The  only 
significant  change  in  this  plot  as  compared  to  figure  6.7  is 
the  movement  of  the  zero  in  the  off  diagonal  pole -zero  plot. 
The      zero  is     seen    to      shift   to      one   of      the  pole      locations 
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Figure  6.9   Open-loop  Transfer  Function  2-1  for  Baseline. 
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Figare   6.10        Transfer  Function  2-1    Optimized  (    3   Var) 
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figure 6.  11    transfer  Function  1-1  Baseline. 
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Figure  6,12   Transfer  Function  1-1  Optimized  (3  variables). 
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Figure  6.13        Closed- loop   Pole-Zero   Plot  Case   3, 
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during  the  optimization.  This  zero  shift  has  the  effect  of 
smoothing  the  frequency  response  curve  in  the  vicinity  of 
the  frequency  of  the  ainimum  singular  value  providing  a  more 
uniform  gain  distribution  around  this  point.  Table  1  shows 
comparison  results  for  the  feedback  gains  for  several  cases 
of   this   tasic   problem. 

The   logical    extension  of   the    problem   to   the   general  case 
is   to   allow     all   four   feedback    gains   to      become    design    vari- 

TA3LE    1 
Comparative  Results   Simple   Problem 


Feedback   Gains 


11 


f.-  variable  f2l=^ 
f    variable  f12=0 
f21  variable  f.-=-5 
All  f's  variable 


0 


1.00006 
0.25254 
.99960 
1.11936 


12 


•51.B3401 

0.0 
50.0 
•60.4718 


21 


0.0 

■0.01071 
0.00025 
0.00280 


22 


0.9998 
1.52332 
0.99038 
0.75478 


ables.  The  pole  placement  and  robustness  routine  can  then 
use  full  freedom  in  choosing  all  four  of  these  feedback 
gains  to  compensate  for  any  cross-coupling  effects  within 
the  system.  Based  on  the  previous  analysis  the  two  diagonal 
gains  would  be  expected  to  approach  1  while  the  upper  off 
diagonal  gain  moves  to  -50  and  the  lower  one  moves  to  0.  The 
flexibility  of  the  ACS  program  was  required  for  the  four 
design  variable  study.  Several  runs  were  made  with  various 
combinations  of  starting  parameters  for  the  feedback  gains 
and   optimizer      routines   before      a   good      design  for      the   case 
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incorporating  all  fcur  feedback  gains  as  variable  was 
obtained.  While  this  points  out  one  of  the  limitations  of 
optimization  design  routines  the  program  was  able  to  develop 
an  improved  design  ever  the  baseline  while  employing  four 
design  variables.  In  the  design  produced  for  this  case  the 
singular  value  level  was  placed  above  0.6  for  a  gain  margin 
of  -4  db  to  9  db  and  phase  margin  of  ±  35  degrees.  In  this 
formulation  the  optimizer  was  able  to  place  the  poles  and 
meet  the  design  singular  value  level.  The  feedback  gains 
produced  by   the    optimizer  were: 

1.11936  -60.4718 

0.0028  0.75478 

which  are  approaching  the  analytic  design  gains.  After 
obtaining  the  feedback  gains  the  OPTSYS  program  was  used  to 
obtain  the  necessary  data  to  'do  a  closed-loop  pole-zerc  map. 
This  plot  is  shown  in  figure  6.14  These  plots  indicate  a 
similar  pole-zero  location  to  that  found  in  the  previous 
three  design  variable  problem.  Again  the  gain  in  the 
affected  channel  has  teen  reduced  to  compensate  for  the  high 
cross-coupling    perturbation    within   the  system. 

The  design  studies  presented  to  this  point  have  been 
based  on  breaking  the  system  loop  at  the  input  as  shewn  in 
figure  6.15.  In  multivariable  theory  the  location  of  the 
break  in  the  loop  changes  the  return  difference  for  the 
system  and  the  transfer  function  formulation.  In  the  figure 
number  1  depicts  a  system  with  an  input  loop  break  point 
while  number  2  is  an  output  loop  break  point  for  output 
return  difference  determination.  The  return  difference  func- 
tion for  the  point  1  is  written  as  I+Fg  while  the  return 
difference  for  point  2  is  I+GF.  The  baseline  system  not  only 
has  low  singular  values  of  the  input  return  difference 
matrix,  the  lowest  being  0.0706,  but  also  has  low  singular 
values  for  the  output  return  difference  matrix  of  point  2  of 
about   the   same    order.      To  demonstrate   the   versatility   of   the 
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Figure   6.14        Pole-Zero  Plot   Simple   Problea    (4   fariables) 
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Figure   6.15       System   Block  Diagram. 

pole  placement  and  icbustness  routine  a  design  was  made  in 
which  the  output  singular  values  were  specified  instead  of 
the  input  values.  The  pole  placement  and  robustness  routine 
produced  a  design  with  gains  that  provided  a  highly  robust 
system  on  the  output  side  with  gain  and  phase  margin  of  -6 
db   to  and  phase    margin  of   ±60    degrees.    The   design   was  not 

robust  on  the  input  side.  Thus,  designing  for  robustness  at 
one  point  in  the  system  does  not  necessarily  give  robustness 
at   all   points   within   the  system. 

Cne  final  case  that  should  be  discussed  is  that  of 
setting  toth  input  and  output  robustness  criteria  at  the 
same  time.  Excellent  results  were  obtained  for  this  case. 
The  design  routine  placed  the  poles  at  -1.97±.009j  with 
feedback   gains    of 

0.85793  -45.91757 

0.00425  0.87137 

The  input  singular  values  were  raised  to  a  level  of  0.74984 
or      -4.5   db     to      12    db      gain   margin      and      ±43   degrees      phase 
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margin.  The  output  values  were  raised  to  above  0.822  which 
corresponds  to  -5  to  15  db  gain  margin  and  ±4  9  degrees  phase 
margin. 

1c  summarize,  it  can  be  stated  that  the  robustness 
problem  for  this  system  exists  in  the  upper  cross-coupling 
channel  (input  two, output  one).  The  lack  of  robustness  can 
be  discovered  in  two  ways.  The  first  method  is  to  plot  the 
open-loop  Bode  plots  cf  each  element  of  the  transfer  matrix 
and  lock  for  extremely  high  gains  and  bandwidths  relative  to 
the  ether  transfer  functions.  The  second  method  examines 
the  singular  values  of  the  return  difference  matrix  for 
magnitude.  Low  singular  values  correspond  to  low  robust- 
ness. The  pole  placement  and  robustness  design  routine  can 
increase  robustness  by  modifying  feedback  gains  to  reduce 
the  effect  of  cross-coupling  within  the  system.  Observing 
the  gain  modification  made  by  the  pole  placement  and  robust- 
ness rcutin€  the  critical  channel  within  the  system  that 
affects  the  robustness  may  be  determined  from  the  Bode 
plots.  The  pole  placement  and  robustness  routine  feedback 
gain  changes  also  cause  zero  shifts  during  the  robustness 
recovery.  The  gain  on  the  open  loop  Bode  plot  for  the 
affected  cress-coupling  channel  is  adjusted  and  the  clcsed- 
loop  zeros,  as  seen  en  the  pole-zero  diagram,  are  shifted. 
This  zero  shift  is  in  a  direction  which  will  combine  with 
system  pcles  to  smooth  the  freguency  response  diagram  in  the 
vicinity   of   the    mininsum   singular    value. 
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VII.    A    HELICOPTER    STABILITY    PROBLEM 

This  chapter  will' deal  with  a  more  practical  application 
of  th€  numerical  optimization  program.  In  this  problem  the 
combined  pole  placement,  robustness  design  procedure  will  be 
applied  to  the  linear  lateral  dynamic  channels  of  a  CH-47 
helicopter.  The  model  is  a  highly  coupled  two-input  two- 
output  system  that  has  been  studied  for  its  basis  robustness 
characteristics  £Ref.  16].  The  usual  procedure  for  design 
of  highly  coupled  systems  is  to  obtain  a  diagonally  dominant 
closed-loop  system.  Ihis  diagonal  system  will  ba  stable  but 
not  robust  to  cross-feed  parameters.  Sandell/et  al,  produced 
three  designs.  Two  cf  these  designs,  while  meeting  basic 
performance  criteria,  had  poor  robustness.  The  third  design 
was  a  relatively  good  design.  The  numerical  optimization 
technigue  developed  in  the  thesis  was  applied  to  the  two 
poor  designs  and  shown  to  provide  substantial  improvement  in 
robustness. 

The  systems  were  designed  to  satisfy  specifications  to 
step  input  response  and  stability  margins  as  stated  in  mili- 
tary specifications.  Specific  design  parameters  for  each  of 
the  three  designs  presented  were  not  available.  The  designs 
were  all  considered  to  meet  the  performance  specification 
criteria  and  stability  margin  requirements.  It  was  shown 
that  two  of  these  designs  were  extremely  sensitive  tc  model 
errors.  The  classic  Nyguist  techniques  did  not  predict  this 
sensitivity.  The  singular  value  analysis  did  indicate  sensi- 
tivity  problems. 

The  model  is  that  of  a  CH-47B  helicopter  lateral  dynamic 
system   in  hover.    The   dynamic    model   of    the    system    is 
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x    ■    A   x    ♦   3    u 


(7.1) 


X    =(v,p,r,y) 


(7.2) 


u 


=   (SS'Sc) 


(7.3) 


where 

2.27 

-1.42 

-0.15 

31.99 

A  = 

0.01 

-0.7 

-0.  07 

0 

0.04 

-0.05 

-0.  5 

0 

0 

1 

0.  11 

0 

and    where 

"  0.  12 

0.95 

B   = 

0.04 
.34 

-8.37 
0.02 

0 

0 

with    full   state    available  for   feedback.    Table   3   is   a   summary 
of   parameters.       The    system    is   not   open-loop   stable. 

Three  control  laws  are  formulated  to  satisfy  the  desired 
perfcriance  specifications.  Equation  7.4  is  the  basic 
contrcl   law. 


u    =   -F  «x   +   h-.   <fc 
i=  1,  2  or   3 

where  the  following  values  of  I   and  h  were  are: 


(7.4) 


f,  = 

' -1.72 

-23.5         70.6 

5  95. 

0.024   - 

-2.71         0.368      ■ 

-7.99 

*l  = 

0.198 

154.0       18.3 

142. 0" 

.    -0.01 

-1.592    -0.  189 

-1.47 

FJ = 

» 
0 

0                25.5 

0 

h  o 

-4                   0 

-27 
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TABLE    2 
CH-46   Helicopter    Parameter  Definitions 


Variable 

Units 

Description 

V 

ft/sec 

Vehicle  body  y-axis  earth  relative 

velocity 

component 

P 

rad/sec 

Roll  rate 

r 

rad/sec 

Yaw  rate 

* 

rad 

Roll  attitude  angle 

6B 

in 

Yaw  rate  rotor  deflection  control 

6c 

in 

Roll  rate  rotor  deflection  control 

■  i        ■      '  -   ■  ■ 

h.  ■ 


hz  = 


hs  = 


The    tfc    is   a  step  input     command   that  must   track.      Figure 

7.1  is  a  diagram  of  the  control  structure  for  the  feedback 
contrcl  laws.  Figure  7.2  is  a  detailed  layout  of  the 
system.  All  three  cf  the  designs  have  negative  real  eigen- 
values and  jrovide  stable  overdamped  responses.  All  three  of 
the   above   designs   meet   the    desired   design    specifications. 

The  loop  Bode  clots,  [Eef.  16],  indicate  all  three 
designs  to  have  good  stability  margins  when  considered  a 
loop   at   a   time.      Since   these      designs   were   all   full   state  LQ 
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Figure  7.1    Feedback  Control  Structure, 
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Figure  7.2    System  Diagraa. 
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designs  Kith  a  diagonal  control  weighting  matrix  they  should 
posses  at  least  -6  dt  to  infinite  gain  margin  an!  60  degrees 
of  phase  margin.  fchen  the  singular  values  of  the  three 
designs  are  computed  a  robustness  problem  is  indicated  by 
low  singular  values  cf  two  of  the  designs.  Figure  7.3  pres- 
ents singular  value  plots  of  all  three  designs.  This  plot 
shows  that  designs  1  and  2  both  have  very  low  minimum 
singular  values  for  the  return  difference  matrix.  Design  1 
goes  as  low  as  -20  db  near  10  rad/sec  in  frequency  while 
design  2  is  down  to  -34  db  at  frequencies  up  to  1000 
rad/sec-  Design  3  is  a  good  design  with  singular  values  that 
remain  above  one  throughout  the  frequency  range  of  interest. 
Using  the  universal  gain  and  phase  diagram  as  discussed  in 
Chapter  6  this  equates  to  a  gain  margin  of  about  -6  db  to 
infinity   and  a    phase    nargin    of   60   degrees   as   expected. 

The  system  stability  in  design  1  may  be  affected  by 
perturbations  occurring  in  the  actuators  as  shown  in  figure 
7.4.  If  the  output  axis  coupling  from  ^  spills  into  £g  in 
the  frequency  range  from  0.5  to  50.0  rad/sec  with  a  magni- 
tude cf  0.12  and  a  phase  of  60  degrees  then  the  system  can 
become  unstable.  This  could  be  caused  by  nonlinear  terms, 
worn    parts,    or   systen  saturation. 

The  second  design  as  shown  in  figure  7.5  may  have  a 
stability   derivative    variation      between  Sq   and  p.  If   this 

derivative  varies  from  about  0.04  to  -0.96  the  system  can 
become  unstable.  The  center  of  gravity  location,  trim  of  the 
aircraft   and  rotor   coupling    can   all   affect   this  ierivative. 

The  pole  placement  and  robustness  design  technique 
developed  for  this  thesis  research  was  applied  to  designs  1 
and  2  tc  obtain  an  improved  robustness  for  these  designs. 
Design  1  will  be  considered  first.  In  this  design  a 
cross-feed  perturbation  through  the  actuator  can  produce 
instability  as  indicated  by  the  low  singular  values  at  lower 
frequencies.      To   study    this    problem      with    the    pole    placement 
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Figure  7.3   Alphatech  Design  Singular  Value  Plots. 

88 


^  (l/2-=GJ<50), 


Figure   7.4        Design    One   Perturbation  In pat. 
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Figure   7.5        resign   Two   Perturbation   Inpint. 

90 


and  robustness  design  method  it  is  assumed  that  the  eigenva- 
lues (pcle  locations)  of  the  system  as  developed  in 
[Ref.    16]-  are  the     required      poles      for   the     performance 

criteria.  Once   the     pole    locations      are   set,        robustness 

criteria  must  be  selected.  From  the  universal  gain  and  phase 
margin  curve  discussed  earlier  two  choices  of  singular  value 
levels  were  made  for  this  problem.  The  first  singular  value 
level  chosen  was  0.6.  This  corresponds  to  a  gain  margin  of 
-4.0  db  to  8  db  and  a  phase  margin  of  about  35  degrees.  The 
second  value  chosen  was  1.0.  With  corresponding  gain  margin 
of  -6  db  to  co  and  a  phase  margin  of  60  degrees  which  are 
the  characteristics  of  a  LQ  regulator  design  with  diagonal 
weighting   matrices. 

For  singular  value  level  0.6  the  pole  placement  and 
robustness  design  routine  places  the  poles  as  shown  in  table 
2  The  slight  differences  in  these  pole  locatioas  appear  to 
have  insignificant  effect  on  the  performance  as  shown  in  the 
response  curves  of  the  system.  The  feedback  gain  adjustment 
moves  the  minimum  singular  value  from  about  0.11  with  very 
poor  phase'  and  gain  margins  to  a  level  of  0.66.  This  is 
above  the  desired  values  of  gain  and  pnase.  The  improvement 
in  robustness  came  from  modifying  the  feedback  gains  in 
channel  £g .  3y  greatly  reducing  the  gains  in  channel  S&  the 
optimizer  minimizes  the  influence  of  the  cross-coupling  from 
channel  Sc  -  In  this  way  a  much  larger  spill  over  of  channel 
Sc,  may  be  tolerated  through  the  actuator  before  the  system 
will    become   unstable. 

The  primary  mechanism  of  robustness  improvement  in  this 
problem  was  a  reduction  in  the  gain  levels  of  the  affected 
channel.  The  feedback  gains  presented  in  table  4  show  the 
modification  of  these  gains  from  those  utilized  in 
[Ref.  16].  While  all  the  gains  are  modified  the  f;^  gain 
undergoes  a  significantly  larger  change  than  the  ether 
design    1    case      1    feedback  gains.         Looking      at   the   open-loop 
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TABLE    3 
Design   One    Pole   Placement 


Pole 

Desired  Location 

Actual  Location 

1 

-24.7977 

-24.7893 

2 

-11.3635 

-12.0083 

3 

-10.3288 

-10.7752 

4 

-  2.1005 

-  2.1181 

transfer  furctions      cf   this   optimized    problem     aa d  comparing 
them   with  a   non-optimized  set      of   transfer   functions   for   0& 
to  as      shown  in      figure   7.6      and   7.7      show  very      little 

change  in  the  system  gain,  both  are  around  42.0  db,  and  a 
bandwidth      of   abcut      20    rad/sec.  There      is   a      significant 

change  in  the  phase  diagram  which  is  caused  by  the  zero 
location   shift. 

The  transfer  function  for  $c  to  £c  depicted  a  gain 
increase  of  about  3  db  for  the  optimized  design  while  the 
phase      remained    similar     for   both      transfer   functions.  In 

transfer  function  Sc  to  Sa  some  important  aspects  of  the 
problem  are  observed.  The  Bode  diagram  of  the  open-loop 
transfer  function  of  &c  to  £&  clearly  indicates  the 
cross-coupling  problems  and  the  pole  placement  and  robust- 
ness design  routine's  mechanism  of  optimizing  the  system 
gains  to  increase  robustness.  The  gain  is  reduced  from  93 
db   to      about  78    db   and      the    bandwidth   is   reduced      from    above 
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Figure   7.6        Transfer  Function  £&-£&  Design  One,    Nonoptiaized. 
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Figure   7.7       Transfer  Function  d&Ss  Design  One  #    Optimized. 
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100  rad/sec  to  about  30  rad/sec.  The  reduction  of  bandwidth 
and  gain  in  the  loop  Sc  to  $&  yields  an  increased  toler- 
ance to  perturbation.  The  transfer  function  for  /c  to 
<$3  is  shown  in  figures  7.8  and  7.9.  In  the  transfer  func- 
tions for  £&  to  £c.  the  bandwidth  is  slightly  increased 
from  0.6  to  0.8  rad/sec  and  the  gain  actually  increased  from 
9      dt   to      14  db.  The  big      change  in      the   overall      system, 

however,  is  in  the  transfer  function  from  £c  to  £a  .  This 
is  the  channel  that  the  destabilizing  perturbation  enters 
and  by  greatly  reducing  the  gain  and  bandwidth  in  this 
channel  through  a  change  in  feedback  gain,  the  optimizer 
routine  has  brought  the  entire  system  gains  to  more  balanced 
conditions    and    recovered  a    highly   robust   design. 

The  gain  changes  associated  with  the  robustness  improve- 
ment cause  the  zeros  ox  the  various  closed-loop  pole-zero 
diagram  tc  move.  A  comparison  of  the  eight  pcle-zero 
diagrams  is  shown  in  figures  7.10  to  7.13.  The  significant 
feature  of  these  pcle-zero  diagrams  is  the  shift  of  the 
zeros  of  the  optimized  design  in  a  direction  that  attempts 
to  egualize  or  balance  the  frequency  response  for  frequen- 
cies in  the  vicinity  of  the  minimum  singular  values.  The 
pole-zero  diagram  of  $b  to  v  will  be  discussed  as  an  example 
of  this  effect.  In  7.10  the  nonoptimized  zeros  are  located 
about  -2  ±  4. 5 j  and  -19.8.  When  the  pole  placement  and 
robustness  routine  has  completed  the  feedback  gain  modifica- 
tion these  zeros  have  shifted  near  -11  and  -6±3j  .  The 
effect  of  these  zero  shifts  is  to  combine  with  the  pcle 
locations  tc  egualize  the  frequency  response  as  depicted  in 
figures  7.14  and  7.15.  Zero  shifts  for  the  remainder  of  the 
transfer  functions  provide  similar  results  in  the  ether 
channels.  By  moving  toward  the  frequencies  associated  with 
the  minimum  singular  values  the  zeros  have  balanced  the 
overall  freguency  response  of  the  system  in  aach  channel. 
While   the   channel   gain   modification   is   the   primary    mechanism 
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Figure  7.8   Transfer  Function  $l~£b   Design  One,  Noncptimized 
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Figure   7.9        Transfer  Function  $c  " S3  Design  One,    Optiaized. 
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Figure   7.10        Pole- Zero  Plots   for   Design    1. 
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Figure  7.11    Pole-Zero  Plots  for  Design  1  (cont.) 
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Figure  7,13    Pole-Zero  Plots  for  Design  1  (cont.) 
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Figure  7.14    The  $g   to  7  Frequency  Response,  Nonopt.. 
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for  robustness  recovery  the  zero  shift  associated  with  the 
feedback  gain  changes  is  directly  related  to  the  overall 
frequency  response    of   the  system. 

In   design    1    case    2   the      minimum  singular    value    level   was 
increased   from    0.6    to   1.0.  The  pole   placement    and   robust- 

ness   design   procedure     also   gave   a   better   design      for  case   2 

T1BLE   4 
Helicopter  Problei  Feedback  Gains 


DESIGN 

GAIN  VALUES 

f   f 

f11  f 

LI21  2 

12  £l3  £l4 

22  x23  J24 

One 
Case  1 

-14.77726 
-0.00567 

2.15858 
-2.556^6 

77.96629 
0.39039 

-32.91595 
-15.04805 

One  • 
Case  2 

-6.36249 
0.00813 

-1.53746 
-2.60073 

72.85013 
0.65395 

74.20387 
-13.47128 

Two 

-0,96660 
0.00916 

-4.54597 
-3.03069 

0.75436 
-0.36391 

11.96808 
-0.34991  ■ 

by  modifying  the  feedback  gains  as  shown  in  table  4.  The 
change  in  the  transfer  functions  for  Sc  to  &g  (figure 
7.16  and  7,17)  is  similar  to  the  change  seen  in  figures7.8 
and  7.9.  The  gain  was  reduced  from  above  90  db  to  about  75 
db  and  the  bandwidth  cut  from  above  100  rad/sec  to  about  25 
rad/sec.  The  added  requirement  of  increased  robustness  did 
not  result  in  a  significant  change  in  performance. 
[Eef.    16]. 
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Eesign  two  which  is  the  design  that  is  subject  to  param- 
eter uncertainty  in  the  eg  to  p  channel  was  also  studied 
using  the  pole  placement  and  robustness  design  routine.  In 
this  design  the  singular  values  were  very  low  until  frequen- 
cies above  100  rad/sec  were  reached.  To  improve  the  design 
it  was  assumed  that  the  pole  locations  that  corresponded  to 
[Ref.  16]  were  the  required  pole  locations  and  that  a 
robustness  singular  value  level  of  0.6  would  provide 
adequate  gain  and  phase  margin  for  the  design.  Ihe  pole 
placement  and  robustness  routine  adjusted  the  gains  in  this 
problem  until  a  minimum  singular  value  of  0.6  was  obtained. 
During  this  adjustment  the  gains  in  channel  $$  were  consid- 
erably reduced  to  offset  the  cross-coupling  between  the  two 
channels.  Again,  in  this  design  as  before  the  channel 
cross-coupling  between  £c.  and  £g  shows  a  marked  change 
in  bandwidth  and  gain  from  above  100  to  35  rad/sec  and  about 
78  db  to  22  db  respectively.  Figures  7.18  and  7.  19  shew  the 
transfer  function  plots  for  this  term.  Figure  7.20  shows 
the  singular  value  improvement.  In  figure  7.21  the  time 
response      is  plotted.  As      can   be     seen      in     the    plot      the 

improvement  in  robustness  for  this  problem  results  in  very 
sluggish  response   and  degraded  performance. 

To  summarize  for  this  problem,  a  given  performance  level 
has  been  chosen  in  terms  of  pole  locations.  The  level  of 
robustness  has  been  set  for  a  desired  gain  and  phase  margin 
based  on  the  universal  gain  and  phase  margin  curve.  The  pele 
placement  and  robustness  routine  has  been  able  to  improve 
the  robustness  level  by  changing  the  feedback  gains  that 
affect      the     channel    £c.    cross-coupling.  This      robustness 

recovery  is  affected  by  modification  of  the  system  feedback 
gains  in  such  a  manner  that  cross  coupling  gains  are  reduced 
so  that  small  cross-ccupling  perturbations  do  not  drive  the 
system  into  instability.  The  open-loop  transfer  function 
plots   have   teen    used    to   indicate    how    this    mechanism   operates 
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Figure  7.20        Singular   Value  Plot  Design  Two. 
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and  have  teen  shown  tc  .be  an  alternative  indicator  of  chan- 
nels that  may  be  affected  by  cross-feed  perturbations.  The 
pole-zero  diagrams  of  the  closed-loop  transfer  functions  of 
the  transfer  matrix  further  indicate  that  zero  movement  is 
in  a  direction  that  equalizes  the  gain  level  of  the 
frequency  response  curves  in  the  vicinity  of  the  lowest 
singular   values    providing  a    more   balanced   system   response. 
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fill.    SIMPLE    OBSERVER 

The   fcle     placement   and    robustness  design     pcocedure   can 
also   be      used   for      robustness   recovery      in    observer      design. 


L*Oir 


X 


X  =  AX  +  BU 


Y=CX 


/s  ^ 


X  =AX  +  KY  +  BU 


Figure  8.1        Simple  Observer. 

Given  an  observer  as  represented  in  figure  8.1  it  has  been 
shown  that  the  system  differential  equation  may  be  written 
as   eguation   8.  1 


U  -BF 

KC    Ac-BF 


I3 


(8.1) 


At    ■    A   -    KC 


where  x  is  the  state  and  x  is  the  observer  variaOle,  and  the 
solution  may  be  separated  to  independent  solutions  for  the 
feedback  gains  and  the  observer  gains.  This  separation 
allows      the      use   of      the     feedback     gains      to  set      the     pole 
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locations  and  then  a  second  optimization  run  with  set  feed- 
hack  gains  so  that  the  observer  gains  may  be  computed  to 
adjust  tie  system  robustness  level. The  observer  pole  loca- 
tion could  also  be  placed  using  the  observer  gains,  K.  In 
the  currently  implementation  of  the  pole  placement  and 
robustness  design  routine  the  observer  poles  are  simply 
restricted    to   areas    cf   the    left   half    plane. 

In  this  chapter  a  simple  stable  observer  system  will  be 
analyzed  based  on  [Bef.  17].  Given  the  system,  equation 
8.2, 


x   = 


0 


1 


35 
-61 


0      1 

-3    -4 

y=  [2    i]i   +»i 

where    EfT)=E(n)=0   and   E  {j(£  )7  (?)  )  «B  K*[  it )  f  (r)  )  =  £  ( t-f) 

the    feedtack  law   of   equation   8.3 


(8.2) 


with 


u    =  -  (5C      10)  x  ♦    50r 


(8.3) 


An  analysis  has  been  done  to  compare  results  of  the  numer- 
ical optimization  procedure  with  results  presented  in 
[Ref.    17].  The  results      for   an      optimal   regulator      design 

using  guadratic  cost  criteria  as  stated  in  equation  8.4  are 
given   in   table   5 


-  /uv 


Hx+    u2)    dt 


(8.4) 


with    H   =    4  ^5    (  -j/35       1) 

Figure  8.2  shows  a  Nyguist  plot  of  the  full  state  regulator, 
the  optimal  filter  and  a  fast  filter.  The  full  state  design 
had  poles  at  s=  -7. 0±j  2.0  and  feedback  gains  of  50  and  10. 
The  optimal  filter  as  shown  in  table  5  had  poles  of  -7. 0  ±  j 
2.0  with  gain  and  phase  margins  of  -6.75  db  and  ±15  degrees. 
The    optimum   filter      gains   were    30   and   -50.  Using   a   faster 
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Filter  with 

Poles  at   

-22+17. 86 j  ;> 


Optimal  Filter 


Full  State 


Figure   8.2         Nyguist   Plot. 

filter   also     gives    poor   gain      and    phase   margins.  The   gain 

margins  are  on  the  order  of  -.98  while  the  phase  margin  is 
less  than  10  degrees.  The  bandwidth  also  increased  from  12 
to  40  rad/sec.  A  recovery  procedure  based  on  a  modification 
to  the  process  noise  matrix  [Ref.  17]  may  be  applied  to  the 
problem.  This  procedure  can  recover  a  large  amount  of  the 
robustness  that  was  lest  with  the  observer  addition.  Figure 
8.3  shows  data  obtained  for  several  trials  of  the  fictitious 
noise      procedure.  The     gain,  phase      margin      and      ether 

parameters    may    be  found   in    table   5. 
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TABLE    5 
Observer    Parameter   Data 


FILTER 
POLES 

GAIN 
HARGIN 

db 

PHASE 

KARGIN 
deq 

ERROR 
COVARIANCE 

e  [(i-iif.-!!)1] 

STATE 
C0VAR1AXCE 
EfwT) 

FUTEH 
CAM 

Optlntl 
LOG  Design 

-7J2J 

-  6.7S 

15 

9.7             -163 
-163            277 

221 
-613 

-613 
2070 

JO 
-60 

Fut  Filter  Adjust- 
■tent  Procedure 

-22J17.9J 

-     .98 

10 

6280       -12200 
-12200        23800 

130 
-613 

-613 

8520 

720 

-HOC 

Fictitious  "olse 
Adjustment  Procedure 

q2  -  100 
q*  •  500 

q'-IO3 

-4.3 

-13.1 

-  7.73 

19 

107           -184 
-184             319 

236 
-613 

-613 
1810 

26. S 

-40.2 

-2.9 
-24 

•10.9 

33 

163           -301 
-301             564 

268 
-613 

-613 
1500 

20.4 
-17.7 

-2.5 
-33 

-13.9 

42 

204           -MS 
-MS            743 

295 
-613 

-613 
1360 

16.7 
-1.1 

-2.1 

-100 

-37 

74 

290           -S70 
-570          1170 

317 

-613 

-613 
1200 

f.9 
M.i 

POPLAR   CASE   I 


POPLAR   CASE   II 


1.00018 
-2.027^2 

13,221*4.8 
32.59^76 


GM=-5.5  dB 
PM=117   deg 

GM=-19   dB 
PM=58   deg 


The  singular  values  of  the  state  feedback  and  optimal 
observer  systems  were  computed  for  comparison  t  ith  results 
produced  using  the  pole  placement  and  robustness  design 
recovery  procedure.  Figure  8.4  indicates  a  loss  in  robust- 
ness represented  by  the  observer  singular  values.  lower 
singular  values  are  less  robust.  For  this  single- input 
single-output  observer  the  Nyguist  diagram  will  be  used  to 
define  gain  and  phase  margins.  Since  a  singular  value  of  1.0 
is  indicative  of  a  linear  quadratic  level  of  gain  ard  phase 
margin, i.e.      GM=-6    db,        and   PM=±60   degrees,    this    was   chosen 
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figure   8.3        Nyguist  Plot   for   Robustness  Recovery. 

as  the  design  level.  The  pole  placement  and  robustness 
routine  was  used  to  recover  robustness  while  setting  the 
pole  placement  at  the  state  feedback  pole  locations  of 
-7±j2.  The  minimum  input  singular  value  level  was  set  at 
1.0.  The  pole  placement  and  robustness  routine  was  also  set 
to  place  the  observer  poles  anywhere  between  -100  and  -2 
that  would  provide  robust  design.  The  plant  poles  were 
placed   at      -7.05   ±j    1.82.      The      feedback    gains  far      this  run 


117 


o 
r 


o 


FULL  STATE 


Bo. 

ZO 

M     * 


O 


.1 


OPTIMAL 

OBSERVER 


l'.O  10.0 

FREQUENCY   -  RAD  /  SEC 


100.0 


Figure  8-4        Singular   falues  of  Observer  System. 
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were  50  and  10.09.  The  optimizer  output  produced  a  filter 
with  a  gains  of  1.00018  and  -2.02742.  The  observer  pcle 
locations  were  -1.99  ±  j  .00  1,  very  near  the  plant  zero 
location.  A  significant  charge  was  produced  in  singular 
values.  Figure  8.5  shows  all  singular  value  curves  plotted 
together.  The  optimizer  solution  for  this  problem  is  well 
above  the  optimal  filter  curve  at  low  frequency.  Figure  8.6 
shows  the  Nyguist  plct  of  the  optimizer  developed  design. 
The  system  has  a  gain  margin  of  -6  db  and  a  phase  margin  of 
117  degrees.  The  most  significant  differences  between  the 
two  designs  being  that  the  observer  poles  are  close  to  the 
plant  zero  locations  and  the  filter  gains  are  much  lower  for 
the    optimizer   solution. 

Using  the      OPTSYS  program      with   the     pole  placement      and 
robustness   routine      computed   gains  as   the      design    parameters 
the    data      for   the  observer      filter   was  computed.         The    error 
covariance    matrix   was  found    to  be 
306.5  -456.5 

-456.5  7C0.1 
These  values  compare  favorably  with  the  trends  established 
in  tatle  5.  The  last  comparison  of  the  pole  placement  and 
robustness  design  recovery  procedure  was  an  analysis  of  the 
time  response  curves.  Figure  8.7  shows  the  comparison  plot. 
The  design  obtained  using  the  pole  placement  and  robustness 
routine   did   not    degrade   system    performance. 

Cne  additional  analysis  was  conducted  that  set  the 
desired  design  parameters  slightly  differently.  In  this  run 
the  pc]e  placement  and  robustness  routine  was  set  to  place 
the  observer  poles  between  -10  and  -100.  The  pole  placement 
and  robustness  routine  was  unable  to  totally  satisfy  this 
requirement.  It  violated  one  of  the  constraints  and  moved 
the  larger  observer  pole  to  -2.25  which  is  near  the  optimum 
pole  location  and  also  the  high  q2  values  of  the  fictitious 
noise    procedure.         The      smaller   pole   was   moved      out    to    -70.8 
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Figure   8,5        Singular   Value  Comparison  Plot. 
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Figure  8.6        Nyguist  for    Computed   Robustness   Recovery. 

which  corresponds  to  the  same  level  of  pole  movement  found 
when  g2  becomes  large.  The  singular  value  levels  were  raised 
significantly  as  shown  in  figure  8.5.  The  state  and  error 
covariance  natrices  were  of  the  same  order  as  the  matrices 
found  in  table  5.  This  design  has  a  phase  margin  cf  58 
degrees.  The  gain  margin  was  about  -19  db.  Even  though  seme 
of  the  constraint  conditions  on  the  design  were  not  met  the 
design   demonstrates    excellent   robustness  recovery . 
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The  pcle  placement  and  robustness  design  routine  has 
produced  a  robust  design  for  this  observer  based  system. 
This  design  is  obtained  by  using  a  numerical  optimization  to 
directly  lanipulate  the  feedback  and  filter  gains. 
Modification  of  the  IQ  functional  equation  as  done  in  the 
fictitious  noise  procedure  is  not  required.  The  pole  place- 
ment and  robustness  routine  solution  for  this  problem  has 
lower  gains  than  those  found  using  the  fictitious  noise 
adjustment.  The  routine  provides  a  good,  direct  methodology 
for  selecting  the  feedback  and  filter  gains  for  a  robust 
observer  design  with  excellent  performance. 
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IX.     BOBOST    OBSEEVER    DESIGN 

This  chapter  will  be  devoted  to  a  short  discussion  of 
the  robustness  recovery  of  a  fourth  order  observer  based 
controller.  The  problem  is  a  helicopter  problem,  [Ref.  6]. 
In  this  case  the  helicopter  model  is  that  of  the  longitu- 
dinal ccntrcl  loop  of  a  CH-47.  The  nominal  model  is  taken  to 
be  the  system  of  equations  9.1  and  9.2  for  an  aircraft  speed 
of   forty  kncts. 


x 
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In  this  problem  the  controller  is  formulated  as  in  figure 
9.1  which  leads  to  an  open-loop  transfer  function  of  the 
form    of   equation   9.3 


K  (s)G(s)  =F(sI-A+EI+KC)-i  KC(sI-A)-iB 


(9.3) 


The  pcle  placement  and  robustness  recovery  procedure  was 
applied   to   this    problem. 

First,  the  standard  full  state  feedback  design  was 
carried  out  using  the  Naval  Postgraduate  versioa  of  OPTSYS. 
This  design  produced  excellent  singular  value  output  for  the 
return      difference    as     shown      in      figure    9.2.  The      lowest 

singular  value  being  essentially  1  ,  corresponding  to  a  LQ 
design  with  -6  db  to  co  gain  margin  and  60  degrees  of  phase 
margin.    The   time  response  of   the   system   was  good    as    shown   in 
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Figure  9-1        Observer   Based  Controller. 

figure  9.3/  reaching  steady-state  in  about  four  seconds  with 
only  a  slight  overshoot.  Since  the  characteristics  for  this 
system  are  acceptable  no  further  design  iteration  was 
carried  out.  The  full  state  feedback  became  the  baseline 
design.  Assuming  that  full  state  feedback  was  not  available 
and  only  two  measurements  could  be  produced  an  observer  was 
developed  to  control  the  two  measured  outputs,  vertical 
velocity  and  pitch  attitude.  Using  the  measurement  matrix  of 
equation  9.2,  OPTSYS  was  used  to  develop  an  optimal  observer 
system  for  this  problem.  The  singular  values  of  the  return 
difference  matrix  and  the  time  response  were  plotted  for 
comparison  to  the  full  state  design.  These  are  shown  in 
figures  9.2  and  9.3  respectively.  Note  that  the  singular 
value  of  the  return  difference  matrix  is  as  low  as  0.16  at  4 
rad/sec.  This  equates  to  a  gain  margin  of  -1  db  to  1  db  and 
less    than      10   degrees  phase      margin.        The   time      response  is 
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Figure   9.2        Singular  Value  Plot   of  Observer  Results. 
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Figure  9-3    Time  Response  for  Observer  System. 
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plotted   in   figure  9.3-      The    steady-state   was   reached   for    the 
velocity   after    a   slight   overshoot   in   about   four   seconds. 

As  the  final  step  in  the  analysis  the  pole  placement  and 
robustness  routine  was  employed  to  recover  robustness  cf  the 
observer  based  system.  The  pole  placement  and  robustness 
routine  was  first  used  to  set  the  poles  at  approximately  the 
same  location  as  the  poles  of  the  LQ  regulator.  These  were 
assumed  to  he  the  desired  pole  locations.  The  pole  placement 
and  robustness  routine  was  then  used  to  vary  the  filter 
gains,  K,  until  the  desired  level  of  robustness  was  reached. 
The  desired  singular  value  level  chosen  was  1.0  which  corre- 
sponds tc  -6  db  to  co  db  gain  margin  and  60  degrees  of  phase 
margin.  Figure  9.2  shows  that  the  pole  placement  and  robust- 
ness procedure  failed  to  totally  recover  the  robustness 
level  to  1.0.  The  minimum  singular  value  reaci ed  was  only 
0.936.  Ihis  singular  value  equates  to  a  gain  margin  of  -5.0 
db  to  co  and  a  phase  margin  of  55  degrees.  While  this  is 
slightly  less  than  the  design  objective,  it  is  Ear  superior 
to  the  optimal  observer  design  discussed  in  the  previous 
paragraphs.  This  recovery  was  made  by  makicg  an  optimiza- 
tion run,  finding  the  filter  gains  with  near  zero  values  and 
freezing  these  values  at  zero.  This  reduced  the  number  of 
design  variables  the  pole  placement  and  robustness  routine 
was  required  to  manipulate  in  a  second  optimization  run  and 
gave    a   higher   robustness  solution. 

This  analysis  and  the  second  order  observer  analysis 
presented  earlier  clearly  indicate  that  it  is  possible  to 
use  the  pcle  placement  and  robustness  procedure  and  separa- 
tion principle  to  develop  a  robustness  recovery  procedure. 
The  pole  placement  aud  robustness  routine  provides  robust- 
ness recovery  by  direct  modification  of  the  feedback  and 
filter  gains.  This  procedure  requires  no  modification  to  LQ 
cost  functicnals  or  ether  parameters  as  done  in  the  ficti- 
tious   noise     adjustment   method   commonly    used     for    robustness 
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recovery.  By  providing  direct  gain  adjustment  the  pole 
placement  and  robustness  procedure  results  in  a  practical 
design  with  relatively  low  observer  gains  and  good  perform- 
ance. The  procedure  is  simple  and  straight  forward  with  the 
only  difficulty  being  a  requirement  to  sometimes  modify 
initial  starting  values  or  optimizer  codes  to  force  the 
solution  toward  the  desired  point. 
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X.    CONCLUSIONS 

An  effective  method  of  robustness  multivaria ble  control 
design  utilizing  a  numerical  optimization  based  algorithm 
has  teen  developed.  The  pole  placement  and  robustness  design 
routine  coupled  with  the  Automated  Design  Synthesis  program 
provides  the  designer  an  excellent  tool  with  which  tc  attack 
the   robust    design   problem. 

The  pole  placement  and  robustness  design  routine  has 
demonstrated  the  capability  of  providing  designs  that  solve 
the  problems  caused  by  cross-coupling  perturbations  which 
reduce  robustness  in  mult ivariable  systems.  This  design 
improvement  is  accomplished  by  modifying  the  system  feedback 
gains  in  such  a  manner  that  the  gain  in  channels  that  are 
affected  by  cross-coupling  perturbations  is  equalized  with 
other  system  gains  to  reduce  this  cross-coupling  effect.  The 
gain  changes  are  accompanied  by  zero  shifts  which  also 
influence  the  gain  distribution  and  frequency  response  of 
the   system. 

Perturbation  problems  in  multivariable  systems  have  been 
shown  to  be  detectable  by  singular  value  analysis  and  by 
using  the  Eode  magnitude  diagram  of  the  open- Id  op  transfer 
functions  of  the  system.  In  the  open  loop  transfer  function 
large  differentials  in  Bode  gains  and  bandwidths  are  indica- 
tive of  problem  areas  for  cross-coupling  perturbations. 
Robustness  is  obtained  by  the  pole  placement  and  robustness 
design  program  by  modifying  those  gains  and  bandwidths  asso- 
ciated with  the  cross-coupling  perturbations  thus  reducing 
the  amount  of  energy  coupled  from  the  perturbation  into 
other  channels.  An  associated  zero  shift  has  been  observed 
when  these  gain  modifications  take  place.  This  zero  shift  is 
in    the   direction      of   poles    that   are  located      in   the   vicinity 
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cf   the    frequency   of    the    minimum      singular    value  and    tends   to 
equalize   the  frequency   response   curve   gains   in  this   region. 

The  use  of  numerical  optimization  to  recover  robustness 
in  observer  based  designs  was  demonstrated.  The  pole  place- 
ment and  robustness  routine  was  applied  to  problems 
previously  solved  using  the  fictitious  noise  procedure  for 
robustness  recovery.  The  direct  manipulation  of  feedback  and 
filter  gains  by  the  pole  placement  and  robustness  routine 
provided  a  highly  robust  design  with  relatively  lew  filter 
gains.  The  problem  cf  robustness  recovery  in  filter- 
observer  designs  has  been  solved  in  a  straight  forward  and 
highly   practical   manner. 
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